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1. SUMMARY 


A new method of analysis for predicting thrust augmenting ejector characteristics is 
presented and the results of computations for a test case discussed. The impetus for the 
development of the method is based on three simulation requirements: (1) a predictive 
analytic procedure is needed, (2) the method must accommodate some form of turbulent 
flow characterization and interchange, and ( 3 ) the final system of equations must be 
amicable to a real-time simulation objective. A literature search revealed an absence of 
ejector research consistent with these combined objectives and that was capable of 
describing transient flows. 

Within the general framework of the control-volume formulation for continuity, 
momentum, and energy there exist time derivatives of field variable volume integrals. Since 
these volume integrals cannot be converted into surface integrals, an approximation for the 
field variable spatial distribution must be made. Under the assumption that the ejector 
mixing region physics dominate ejector performance characteristics, spatial sub-division of 
the mixing region permits each sub-volume to be approximated by characteristic velocity, 
pressure, ana temperature profiles. A description of turbulent flow is provided with 
Abramovich-type self similar turbulent flow field variable profiles. Time derivatives of the 
volume integral reduce to time derivatives of the field variable characteristics, and, with 
treatment of the surface integrals in the "usual" way, a set of differential equations in time 
evolves. With the intent to focus primarily on results for ejector thrust, very few (less than 
ten) subdivisions of the mixing region are needed, and, therefore, a terse description of the 
ejector mixing region is obtained. Although a crude description of the turbulent ejector jet 
interaction is employed, the final system of equations can potentially provide real-time 
thrust predictions, thereby meeting the aforementioned objectives. 

Since a step-change in the ejector driving nozzle flow is representative of typical 
ejector operation, an example prediction of this situation is employed as a test case for 
application of the methodology. 



2. EJECTOR SIMULATION PERSPECTIVE 


Role of the Simulation 

Research on design methodologies for integrated aircraft and propulsion flight 
control systems requires accurate subsystem component simulations. In principle, these 
simulations must mimic steady-state and transient component effects. A NASA Lewis 
research program is currently underway to develop a "real-time" simulation - including 
system transients - for Short Take-Off Vertical Landing (STOVL) aircraft. Thrust 
augmenting ejectors are considered potentially valuable propulsion subsystem elements 
for the powered-lift aspect of STOVL aircraft. To explore ejector concepts further, the 
initial NASA Lewis STOVL system simulation requires a thrust augmenting ejector 
sub-system simulation. Unfortunately, an ejector simulation that includes ejector 
transients is not currently available; the purpose of the present work is to develop one. 

An ejector is a mechanically simple fluidic pump composed essentially of two 
components: (1) a "primary" jet nozzle issuing into (2) a shroud. This arrangement 
permits entrainment and acceleration of a secondary flow (within the shroud) by the 
primary jet. A diffuser section attached to the shroud allows control over the ejector 
discharge pressure. Figure 1 illustrates a generic thrust augmenting ejector. 

From a system simulation point of view, the ejector participates in the description 
of the aircraft "plant dynamics" as shown in figure 2. 

Three Competing Requirements 

There are three basic requirements the ejector simulation must meet: 

1. The mathematical model must be predictive (not parametric) in nature. 

2. Some approximation of turbulent flow characteristics inside the ejector mixing 
region must be made. 

3. The final system of equations describing the ejector must be amicable to the 
NASA/Lewis real-time simulation objective. 

To meet the first requirement, only two data sets should be prescribed: (a) the primary 
jet "control-valve" setting (with associated thermodynamic data), and (b) free-stream 
atmospheric properties. From this, the secondary inlet condition and the ejector thrust 
augmentation (as a function of time) are predicted. A parametric method would specify 
all primary and secondary conditions (in contrast to the predictive method where the 
latter is unknown). 

Central to the idea of an accurate ejector simulation, the second requirement 
points out the need for some type of characterization of turbulent mixing and 
entrainment phenomena in the mixing region; a description of these effects is essential 
for transient- type analyses. 

The third requirement emphasizes the ejector simulation must not become a CPU 
bottleneck in the final STOVL simulation. 
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At this point in time these simulation goals are at odds with each other; methods 
relying on first principles for the predictive facet of analysis are accurate, but not 
executable in real-time. 

Perspective on Previous Work 

Ejector research goes back about half a century and has resulted in a multitude of 
papers on various aspects of ejector performance, optimization, and analytic methods. 
Frequently referenced in recent literature is the work of Porter and Squires[1981] 
whose survey produced a compilation of over 1600 research papers on ejectors; much 
more research on ejectors has been done since that review. No attempt is made in the 
present work to extensively discuss the history of ejector research. Rather, it is more of 
interest to note that in the survey conducted for this work (see References) unsteady 
flow research focuses on the pulsed ejector flow problem. Such flows are of interest, for 
example, in chemical laser applications (see Anderson [1970a, 1970b, 1976], Johnson 
[1966], and Petrie, Addy, and Dutton[1985]). No papers were identified that dealt with 
transient thrust augmenting ejector flows in a way consistent with a real-time simulation 
goal. 


Methods of analysis are extracted from either a control-volume or material-vol- 
ume formulation. In the former, a set of algebraic equations is (often) based on 
quasi-one-dimensional and isentropic conditions, while in the latter a multidimensional 
non-linear system of differential equations is obtained. The first generally gives way to 
parametric studies while the second is more predictive in nature. 

Although bases of analysis can be polarized as described above, there are many 
methods of implementation. For any specific perspective on formulation involving, for 
instance, a system of differential equations, details of analysis depend on the analyst’s 
selection of a solution method, e.g., finite-difference, finite-element, or method of 
characteristics. We summarize the extremes of ejector simulation as follows: 

1. A control volume analysis is by far the most straightforward and widely reported 

method, but, since the system of equations is under-determined (when only the 
ambient and primary nozzle conditions are given), some key parameters must be 
specified. Although the system has potential for real-time simulation, the need to 
prescribe, for instance, mass entrainment ratio yields that (in principle) this is at 
nest a parametric method of analysis. 

2. A material-volume analysis draws on the full Navier-Stokes and energy equations, 

and provides an opportunity for a fundamental introduction of tubulent flow 
characteristics. It is therefore possible to be predictive in nature, but only at a 
considerable CPU expense; these machine computations are not likely to fall 
within a real-time framework. 

Present Approach 

In order to generate a predictive method of analysis an empirically based model 
for the turbulent interaction region is explored within the framework of a control 
volume analysis. This approach provides a rational foundation for the introduction of 
steady flow data to "calibrate" an unsteady flow simulation; there is no intent here to 
provide a multi-dimensional CFD code based on first principles. 
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Elements of Current Work 

A methodology for simulating thrust augmenting ejector performance is described 
in the present work. Section 3 describes some typical ejector approximations and 
considerations related to ejector operation and operating regimes. Since the final 
system of equations are extracted from a control-volume formulation, control-volume 
equations for an arbitrary ejector control volume are given in Section 4; this enables the 
general (time-dependent) surface velocity to be correctly introduced into the system. 
Section 4 also expands on the simplification of one-dimensional flow and application of 
the equations to the inlet, mixing region, and diffuser. 

Section 5 remarks on the simplification offered by a steady-state analysis and 
provides a descriptive solution procedure and results from the same. Section 6 looks at 
three approaches for unsteady flow analysis and details the Finite- Volume adaptation 
used in the present work. 

A test case is examined in Section 6; this leads the the concluding remarks in 
Section 8. 

Appendicies A-J provide extensive detail on points of analyses traditionally 
assumed "intuitively obvious". Such information is contained in the present work for (a) 
completeness of documentation of the proposed method, (b) the capability to repro- 
duce the derived results, and (c) illustration of inadvertently implicit assumptions. 
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3. PRELIMINARY CONSIDERATIONS 


Qualitative Ejector Characteristics 

A thrust augmenting ejector is often described as a fluidic pump that employs the 
momentum of a nigh velocity jet from the primary flow nozzle (drive flow) to entrain 
and pressurize a secondary (suction) stream; a typical thrust augmenting ejector consists 
of four basic components: 

1. a high pressure nozzle to accelerate the primary flow, 

2. an inlet section to accelerate the secondary flow 

3. an intermediate mixing section to permit momentum exchange between the 
primary and secondary flows, and 

4. a diffuser to match the discharge pressure (static) with the ambient. 

Overviews of the characteristics for this general ejector configuration have been 
given recently by (among others) Koenig et. al.[1981], Minardi [1982], and Bevilaqua 
[1984], Also, the proceedings of the 1981 Ejector Workshop for Aerospace Applications 
(Braden et. al. [1982]) covers many issues in ejector technology and simulation. It is well 
known that the (irreversible) mixing of the primary and secondary streams results in a 
local static pressure that is less than the ambient; this is the origin of the suction effect 
on the secondary stream. Recovery of the stat ic pressure in the diffuser results in a net 
thrust component from the difference of the (integrated) pressure distribution of the 
diffuser and inlet 1 . It is therefore of interest in ejector design to contour the inlet and 
diffuser so as to maximize the suction effect and diffuser pressure recovery. A typical 
ejector wall pressure distribution is shown in Figure 3 (taken from Bernal and Sarohia 
[1983]; also see the work of Miller and Comings [1958]). Minimization of nozzle drag 
(direct and ram) is also important, as is the need to avoid shroud leakages. An 
approximate relation for the ejector system thrust is given by 
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m 


\ pit \ p 


+ p(I • n)d A 
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where the first term is an approximation for the primary nozzle thrust, the second the 
net surface pressure integral, and the last the sum of viscous and pressure losses. The 
obvious design goal is to have an ejector where the third term is minimized. In analysis, 
the object is to predict the velocities and pressures such that all terms in the thrust 
equation can be evaluated; as a practical matter this is not an easy task and is the 
impetus for the variety of approximate methods that exist. For example, in a 
steady-state analysis the total thrust of the system can be computed from application of 
conservation of momentum and mass for a control volume corresponding to the duct 
boundary - in this case there is no need to integrate the pressure distribution over the 


1 Analogous to the theory of lift on a wing in an inviscid flow, the inviscid thrust in an ejector is attributed 
to the net circulation that arises when the flow streamlines in the shroud arc directed longitudinally from 
what would otherwise (in the absence of the shroud) be at an angle to the centerline. 
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wall of the shroud. At the other extreme finite-element or finite-difference methods of 
analysis are used to compute detailed flow velocities and pressures along the wall 
bounding the flow. 

It is evident that (as in many other fluid flow problems) ejector analyses suffer 
from a lack of understanding of turbulent flow. Current mathematical descriptions of 
turbulent flow yield non-linear, time-dependent equations; non-unique solutions are 
also an important consideration. Phenomenological descriptions prevail, however, and 
"calibrated" versions seem to simulate flows well. Several ejector studies blend 
empirical results, aero-thermodynamics, and control volume (or numerical) approaches 
quite successfully in the analysis of ejector performance (see for instance, Salter [1975], 
or Tavella and Roberts [1984]). Nonetheless, experimental work continues to improve 
understanding (again, an example is the work of Bernal and Sarohia [1983]). Some 
conclusions from experimental studies assist in the characterization of ejector behavior: 

1. The level of thrust augmentation does not vary noticeably with primary pressure 
ratio (the rate at which the primary jet spreads is not a function of its initial 
velocity). 

2. Thrust augmentation increases nearly linearly with diffuser area ratio, up to an 
area ratio in the vicinity of 1.5. Augmentation levels off and/or decreases as the 
area ratio is increased further - the duct wall half-angle seems to influence the 
particular trend. 

3. Ejector performance is very sensitive to inlet losses and the thrust efficiency of the 
nozzle, although skin friction losses appear to be very small. The extent of mixing 
of the flow discharge (flow skewness) has a significant impact on ejector 
performance. See Belivaqua [1974], 

4. Velocity profiles in the mixing region tend to be self-similar. 

5. Additional thrust augmentation can be realized for installed ejectors with the use 
of end plates - these cause an otherwise three-dimensional flow to be two-dimen- 
sional (the end plates block flow into the separated region and create a drop in 
pressure at the ejector exit, improving performance). 

Reflecting on the last remark, the present work does not focus extensively on some of 
the aerodynamic installation problems of ejectors described, for instance, by Knott and 
Cudy [1986] or Lund, Tavella, and Roberts [1986], though this is not to say that such 
effects are unimportant or can be overlooked in future ejector examinations. Rather, 
two- vs three- dimensional effects are significant and have considerable bearing on the 
interpretation of experimental data. 

Of principle interest here is that secondary mass flow entrainment and the extent 
of mixing (of the two streams at the ejector discharge station) are closely related to 
thrust augmentation, and influenced by (1) ejector geometry, (2) primary jet character- 
istics, and (3) physical fluid properties. 
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Efficiency and Thmst Augmentation Ratio 

Ejector performance is often quantified by computation of the thrust augmen- 
tation ratio, <t>, and the ejector efficiency, 7. As would be expected, an increase in 
thrust augmentation occurs at the cost of a decreased efficiency, so the "optimum" 
ejector balances the two (in accordance with the prescribed mission). In the 
definition of 0 it is useful to use as a reference the isentropic thrust obtained from 
the expansion of the primary jet to the ambient 

. E SYSTEM 

= p (?) 

r PRIM ARY. IDEAL 


Note that when an isentropic reference thrust is used it is easier to compare the 
performance of ejectors with different nozzles. 

Several equations to characterize "efficiency" have been used in the literature: 

1. The nozzle efficiency is given as the ratio of the system thrust and the sum of 
the thrust for the primary and secondary streams under ideal conditions, 

2. The ejector efficiency can be measured as the ratio of the kinetic energy of the 
ejector efflux to the input energy of the primary nozzle, 

3. A ratio of the input momentum to the discharge momentum could be used (for 
the nozzle or ejector), 

4. Base the ejector efficiency on the concept of thermodynamic availability ( see 
Minardi [1982]), 

5. Compute the ratio of the enthalpy change of the mixed ejector flow to the 
ideal change in enthalpy of the primary flow (again, see Minardi [1982] for details 
and several versions on this). 

Item 2 is chosen for the present work and described in more detail in Section 4. 
Characteristic Surfaces 

Mass flow characteristics can be expressed in a general way as a function of 
stagnation pressure and back-pressure ratios as 
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where it is assumed that ejector geometry and fluid properties are known. Similarly, 
the secondary stream inlet Mach number can be expressed as 
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where Pis is the static pressure of the secondary flow at the point of confluence of 
the two streams. 

A significant contibution to ejector analysis is the work of Addy, Dutton, 
Mikkelson and co-workers [1974, 1981, 1986] who provide a clear view of overall 
ejector characteristics through the presentation of three-dimensional surfaces; the 
ordinate and abscissa use the arguments of the general relations given above. An 
overview of the more detailed discussion by Addy, Dutton, and Mikkelsen [1981] 
follows. 

Three regimes of flow can be described with three-dimensional surfaces that 
have the parameters of equations (2) and (3) as axes; these surfaces are shown in 
figure 4. An important feature of the surface is the "break-off curve that divides the 
"supersonic" and "saturated-supersonic" regimes and the mixed flow regime. In the 
latter, <t> is a function of the ambient pressure level while the former are not. 
Although it is assumed the primary nozzle flow is choked, the distinction of the 
breakoff curve is to mark the development of sonic conditions in the mixing region. 
Under certain conditions the secondary stream velocity can reach sonic conditions at 
the inlet, in which case the flow is choked and the mixing region described by a 
saturated-supersonic regime (see figure 4). Subsonic flows are depicted by the mixed 
regime. 

Thrust augmenting ejectors entrain a second fluid at ambient conditions and 
discharge to the (same) ambient a mixed primary and secondary fluid. In this 
situation 

P 15\0 _ Palm (g) 

P ) P.O P 1 P.O 

so a subset of the three-dimensional surface of figure 4 which is of interest is the 
two-dimensional slice shown in figure 5. Because the three-dimensional surface has 
been drawn in a general way, no specific intent to exclude or include the 
saturated-supersonic region has been made in the description of figure 5. The 
secondary flow characteristic surface of figure 6 illustrates that sonic inlet conditions 
are associated with the saturated-supersonic regime. 

Obviously the problem before us is to establish specific numbers for the axes 
of the characteristic curves (as functions of time). Recognize that those numbers 
which are presently available (and supported by experimental data) are generally for 
steady-state ejector operation and therefore provide little with which to remark on 
transient behavior. 
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Remarks on the Analytic Approach 

In the interest of (eventually) realizing a real-time simulation capability, it is rational to 
begin the ejector analysis with the development of (1-D, integral) control-volume 
equations for a partitioned ejector; nomenclature for an arbitrary control volume is 
given in figure 7. This type of analysis forms a general mathematical structure within 
which methods for providing a transient capability and imitation of the turbulent 
interactions can be explored. 

To be sure, there are much more capable frameworks of analysis, but they have 
only been explored for steady flows and example calculations indicate the approaches 
are extremely CPU intensive. Shen et. al.[1981] investigate the high secondary mass 
flow scenario using finite element analysis. Hedges and Hill [1974] discuss finite-differ- 
ence solutions and review the integral boundary layer analyses of significance. 

Considerable attention is paid in Appendix A and B to the development of the 
integral equations since (a) the nature of the time-dependent terms is important, (b) a 
proper account of a control volume moving in space is required, and (c) an open 
discussion of "intuitively obvious" quantities (presented without derivation) is without 
rigor. The development of these equations is predicated on a point of view appropriate 
for ejectors - "appropriate" is given here by those arguments often presented and 
weathered scrutiny in the literature. Some elements of that point of view are discussed 
below: 


Constant Area Mixing Section 

The ejector shown in figure 1 reflects the assumption that the mixing region 
has a constant cross-sectional area; this selection is intentional and made on a 
theoretical basis. An alternate configuration is the "constant-pressure" ejector in 
which it is assumed the area of the mixing region varies such that, in the 
one-dimensional case, the integrated static pressure is constant over the mixing 
region cross-section. For this configuration, however, determining the necessary 
muring section area distribution and the high-probability of off-design operation are 
problems. 

As it turns out for the traditional steady-state control volume analyses, the 
constant-area and constant-pressure formulations predict comparable performance 
(see Dutton et. al.[1982], Dutton and Carroll [1986]). In fact, the constant area 
formulation leads to a doubled value solution in which the mixed flow is either 
subsonic or supersonic; the solutions for the Mach number are related by the 
relation for the Mach number across a normal shock. Selection of a solution is based 
on compliance with the second law of thermodynamics. 

Minardi [1982] provides a discussion of considerations for which the con- 
stant-area geometry is a necessary, but not sufficient, condition for the analysis. This 
leads to the conclusion that all possible solutions for the mixing region will be 
obtained with the constant-area case, even though certain solutions are not likely to 
occur for typical ejector configurations. 
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Adiabatic Ejector Walls 

It is almost a universal assumption in theoretical ejector studies that the 
ejector walls are insulated; this allows the surface heat flux term in the energy 
equation to be dropped. As will be discussed in detail later, this is not immediately 
an assumption that the flow is isentropic. If we consider the arguments on the 
creation of entropy in Appendix B, 
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then entropy changes exist as long as the fluid is subject to viscous dissipation. The 
combination of adiabatic ejector walls and inviscid flow provides: 


from which a convenient steady flow integral arises: 



0 


( 8 ) 


This relation is most often employed in the analysis of the inlet and diffuser. Since 
Bernoulli’s equation can now be used for the description of pressure across those 
regions an alternate use of the entropy balance is realized; generally speaking, the 
entropy balance becomes a condition on the simultaneous solution for the equations 
of mass, momentum, and energy (mechanical). 

Inviscid Interaction Region 

For ejector flow conditions where the static pressure of the secondary flow is 
less than that of the primary flow (Pgj < Ppj), the primary flow expands and 
interacts with the secondary flow to provide the interface boundary shown in figure 
8. Note the formation of an "aerodynamic throat", that is, a minimum cross-sectional 
area downstream of the inlet associated with the boundary of the inviscid interaction 
of the two streams. If the secondary flow is subsonic upstream of this station, then 
the peak flow velocity at the constriction will be at or below sonic conditions. The 
"supersonic" operating regime is described when M$2=l, and the "mixed" region 
given where At very low secondary flow rates the secondary flow is 

effectively "sealed" off from the primary flow. A number of investigators (Anderson 
[1974a, b], Addy and co-workers [1974, 1981]) have taken advantage of these 
"partitioned" flows, as originally characterized by Fabri[1958]. Of interest here are 
the admission of the following assumptions for this "inviscid interaction" region: 

1. The primary and secondary streams are distinct and do not mix between the 

point of confluence and station 2. 
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2. Viscous interactions occur along the interface of the two streams (mechanism 
of energy exchange between the primary and secondary flow). 

3. Each stream is treated as irrotational. 

4. Although the average pressure of the streams is (potentially) different at each 
streamwise station, the local static pressure is equal at the boundary. 

In the work of Addy et. al. [1981] both streams are treated with an isentropic flow 
assumption (mixed and supersonic regimes). Anderson [1974a] allows that the 
secondary flow is isentropic, but determines the primary flow field from the method 
of characteristics (in his treatment of supersonic flow only). Quain et. al.[1984] relax 
the assumption that the secondary flow is isentropic, but the discussion of results 
from that work is not possible at this time (the paper is written in Chinese, abstract 
in English). 

Implementation of an inviscid inte ration region may simplify the analysis of 
the mixing region, but this does not relax the requirement of an iterative solution 
methodology. Nonetheless, the use of this modelling approach offeres some analytic 
simplification and is worthwhile to pursue if the conditions which precipitate the 
approximation can be identified in advance. 

The low-secondary-flow rate regime is only properly treated with a two-dimen- 
sional flow field analysis. As such, the euphoria surrounding a one-dimensional 
approach is mitigated by these types of flows and asks that the computational 
burden be assumed (for the correct analysis) or the degredation in solution accuracy 
announced. 

Thermodynamic Considerations 

Some elementary thermodynamic assumptions are often made in ejector 
analyses; identification of them here is appropriate: 

1. The state principle of thermodynamics provides that, for a gas, a given 
thermodynamic variable can be described in terms of two thermodynamic 
properties. 

2. Below a gas temperature of (approximately) 600°K, air has essentially a 
constant specific heat, and therefore qualifies as a perfect gas; for generality, 
though, an ideal gas would be used 

3. For a perfect gas the entropy anti enthalpy of the gas are expressible as a 
function of temperature only. 

4. Application of Dalton’s law of partial pressure to a constant volume mixing 
process yields the properties of the mixed flow at station 3 (see Appendix F). 
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Rarely does it come to pass the need to distinguish between an ideal and a perfect 
gas, both are described by the equation of state, but the specific heat of an perfect 
gas is constant, but in an ideal gas the specific heat is a function of temperature. This 
minor point can have significant impact on computational speed if temperature is 
unknown and the specific heat then extracted by an iterative method. 

Summary 

In the present simulation effort, an ejector with a constant area mixing region and 
adiabatic walls will be used to provide thrust augmentation through the mixing of two 
ideal gas streams. A phenomenological approach is taken to the turbulent interaction 
problem, consistent with the objective to formulate a system description amicable to 
real-time simulation. 
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4. CONTROL-VOLUME EQUATIONS FOR ONE-DIMENSIONAL FLOW 
Overview 

A control volume description of the ejector is perceived to be the most efficient 
method of analysis for the real-time simulation goal of the present work. This section 
outlines the application of the general control volume equations (developed in 
Appendix A) to the ejector shown in figures 1 and 10. Note the ejector is partitioned 
into an inlet, mixing, and diffuser region. Development of the mass, momentum, and 
energy equations for each component intends to provide three results: (a) a summary of 
the integral form of the equation, (b) a form useful for the transient ejector analysis, 
and (c) the version obtained if a piecewise-constant velocity is assumed at stations of 
inflow and efflux. Prior to those developments some remarks on the one-dimensional 
flow approximation are made below. 

One-dimensional flow approximation 

Compressible flow in channels is often treated in (practical) control volume 
analyses under the assumption that a quasi-one-dimensional flow exists; the basic 
simplification offered is that velocity gradients can occur along (not across) the 
streamwise axis of the channel (longitudinal axis). See figure 9. 

A real-flow velocity profile is not (necessarily) symmetric and would reflect the 
presence of any viscous, blockage, and Reynolds number effects. In the absence of 
separation the real flow in the vicinity of the boundary must be parallel to the wall. In 
the present work these effects can only be accounted for as far as the continuity 
equation allows. 

Analysis of "traditional" one-dimensional flow yields that the accuracy of this 
assumption depends on the axial gradient of the cross-sectional area. If dA/dx is small 
the assumption is well received (For the analysis of a stream-tube the "approximation" is 
exact). It is worth noting the comparison given by Thompson [1972] between nozzle 
data and the result for a one-dimensional flow analysis. 

Particular attention should be paid to the generosity of the one-dimensional flow 
assumption for the expansion of the primary jet issuing from the nozzle. There is 
considerable breakdown in the approximation - and therefore in the accuracy of the 
analysis - when the secondary static pressure is significantly less than that of the primary 
jet. Recognize, of course, the extreme subsonic case (discussed previously) where the 
secondary flow is effectively "sealed-off and the flow field must be treated as 
two-dimensional (recall figure 8). 

Often it is assumed that the one-dimensional flow is also piecewise constant; this 
allows for multiple streams of constant one-dimensional velocity to be considered and 
provides for very convenient forms of the control volume equations. 

Some accommodation of real-flow effects can be provided through a "skewness 
factor" which mitigates the asumption that the discharge velocity from a duct is 
completely mixed. As a practical matter, such an account is particularly useful for 
examination of ejector configurations where experimental data is available. 
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Skewness Factor 

In order to include (partially, anyway) non-ideal mixing characteristics in the 
ejector analysis, there is often introduced a skewness factor for the flow at the exit of 
a mixing region and for the diffuser exit. It is assumed that knowledge of friction 
losses in the shroud and the value of the skewness factor are adequate for the 
characterization of the net effects of non-ideal mixing (even though local flow 
details cannot be extracted). 

Consider the spatially averaged velocity given by 


( 1 ) 


[< v > 7 A) ( 2 ) 


<v> = — vdA 

A, 


if' 


for the definition 


■ (/■' 


dA / 


If v is uniform, the skewness factor is unity; in a non-uniform flow /? > 1. To get a 
feeling for the magnitude of p, Bevilaqua [1974] notes that a typical ejector inlet 
region has a skewness factor on the order of 1 .8, while the skewness factor at the exit 
of a (well-designed) mixing duct is approximately 1.02. 

The skewness factor is constant for self-similar velocity profiles. 

Salter [1975] discusses a theoretical approach to the prediction of the skewness 
factor, based on the turbulent jet theory of Abramovich [1963]. 

Continuity 

For a control volume moving through space with velocity U(t) the continuity 
equation is 


d_ 

dt 



p ( u - {/ ) • ndA 


0 


(3) 


Define at any instant in time the mass contained in the control volume by 



(4) 


and introduce the relative velocity 

v = a-U = vi 


(5) 
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where the unit normal is defined as positive outward from the surface of the control 
volume. Also, define the mass flowrate past station k as 

rh = “ f p At (v- n)dA ( 6 ) 


Emphasis is placed on (a) the use of the relative velocity, and (b) the presence of the 
negative sign. The impact of the latter is that an inflow is described by a positive mass 
flow; since 


(wa)k = u k (-i-i) = ~v k 


cn 


then for an inflow 


rti i 


- / 


pv k dA 


( 8 ) 


and for an outflow 


(«•«)*♦! = V k + i (t ' 0 = U*.| 


(9) 


m 


k * 1 



p v k d A 


( 10 ) 


Returning to the continuity equation, substitution yields 

= ( 1 1 ) 

This has the intuitively desireable result that when the inflow exceeds the outflow, the 
accumulation of mass is a positive quantity. Although superficially this may appear a 
trivial result, reflect on a formulation involving u instead of v in the definition of the 
mass flowrate; the impact on the latter is significant. 

Introduction of the skewness factor requires the assumption of a constant density 
over the cross-section, 

rh = pv k dA = p v k dA = p k <v> A (12) 

J A k J A k 

Integrals dissapear on the right hand side of the continuity equation if the density 
is considered piecewise constant; the mass flux accross station k is 


< is > 



(pvA) k 


(13) 


rfi 


k 


and therefore 

dm, 


dt 


= (pvA) k - ( pvA) kt 


(14) 


Momentum 

Application of the integral momentum equation to an ejector with adiabatic walls 
and where gravitational effects are neglected provides 


— pud V 

dtJv 


pu(u - U) ■ ridA + 

A ~~ ’ J A 


j> n • S' d A - j> n •( p l')d A 


( 15) 


Introduce the relative velocity, then 


dt 


f P(v 

J V 


+ U)dV = 


p ( u -*•{/) u • n d A 


£ a 

£ a 


(p/')dA 


+ <P n - S' d A 

1 A 


( 16 ) 


Define the net force acting on interior surfaces (thrust) to be characterized by the sum 
of the surface integral of the deviatoric stress tensor and the interior pressure forces 


- £n 


S'dA 


~ jn-CpfJdA 


(17) 


It is important to note that no shear stresses are assumed to be acting over the inlet or 
outlet regions. Also note that the area r represents the (fixed) interior surfaces of the 
ejector unit - those are the areas responsible for the exchange of forces (between the 
fluid and ejector surfaces) to provide the system thrust. 

Expand the time derivative in the momentum equation 

— | p(v+U)dV = d f pvdV + d f pUdV (18) 

dtJ v ~ dtJv dtJv 

and recognize 
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jtL p ^ dv - 


( 19 ) 


Recalling that the reference velocity U is constant over the surface (but remains time 
dependent), then 


d 

dt 


U.m C J 


U 


dm, 


dt 


m. 


dU 

dt 


( 20 ) 


Also 


- f p(y + L/)u • nd A = - <i> pv(v • n)d A - U <f pv ■ rid A 
J A J A J A 


dm, 

py(y- n)d A + U - 
a dt 


(21 ) 


The momentum equation becomes 


d 

dt 


X 


pvdV 


+ 


rn J 


dll 

dt 



pv(i • n)d A 


x 


ri- ( p / )d A + /• 


( 22 ) 


where A implies those surfaces at the inlet and outlet of the control volume, not the 
internal surfaces used in the definition of F . 


Important Note: The volume integral represents the time rate-of-change of the 
momentum of all particles inside the control volume at an instant in time. It cannot be 
converted into a surface integral and requires some estimate of field conditions in the 
ejector mixing region. 

Some simplification of the momentum equation results from a piecewise constant 
velocity, density, and pressure assumption. In this case 


j> pvfy 


n)dA 


- j> n ■ ( p/')dA 


* 

- Z 1 kC l k d k 

k 


( 23 ) 


(2 1 ) 


and therefore 



pyd V 


m , 


dU 

dt 


Z v -k rtl k 


k ri k A k 

k 


F 


( 25 ) 
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An alternate modification involves the skewness factor. For a uniform cross-sec- 
tional area, 


/ 


pv 2 dA 


pp <v> 2 A 


= p rh <v> 


( 26 ) 


then substitution in the appropriate locations in the momentum analysis above provides 
the desired result. 

Energy 

Ejector analyses typically assume all surfaces are adiabatic, but (as mentioned 
earlier) this is not immediately an implication that any of the transport processes are 
isentropic. With this in mind the control volume formulation of the energy equation is: 


~ f (pH-P)dV 
dt J v 


pU-ndA - <fi pll(a- U)- ndA 
A J A 


( 27 ) 


Note the absolute velocity participates in the definition of the stagnation enthalpy 


u 

H = h + — 

2 


( 28 ) 


From this the ideal gas assumption yields 

pH = pc p T 0 


( 29 ) 


then 

^{pc p T 0 -P)dV = j pU-ndA - £ pc p T 0 y - ndA (30) 

which is the desired integral form. 

As before, a piecewise constant velocity, density, and temperature simplifies the 
equation 

dtfv^ PCpT °~ = n k ) + ( 31 ) 


but does not eliminate the volume integral on the RHS. 
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Entropy 

In integral form the entropy equation acts as a condition on the solution of the 
mass, momentum, and enegy equations; from Appendix A, 


-f 

dt J v 


psdV 


£ 


psv • ndA 


£ 


— q-ndA 
aT ' - 


> 0 


(32) 


Acknowledging that the last term vanishes by virtue of the adiabatic ejector assumption, 
a non-zero entropy balance arises from viscous dissipation. 

If the mass flowrate is introduced 


d_ 

dt 


L 


psdV 


£ 


sd th > 0 


(33) 


then the uniform one-dimensional flow assumption yields 


-f 

dt J v 


psdV 


Y s th ( i ,, • n ) > 0 

k 


(34) 


Summary of Basic Equations 

Four basic control volume equations form the foundation of the analysis of the 
transient ejector: 

Mass 


— {m cv } = Yrn k 
dt ' cvl V 

Momentum 

d f dU f 

-- \ p v vdV + rn cv = -<Ppv(vn)dA 
dt J dt Ja 

Energy 

— f [pc T 0 - P)dV = ipU-ndA - f pc p T 0 v ■ rid A 

dt J v J J a 

Entropy 


- £ n 


(pI)dA + r 
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— f psdV 

dtJv 

For the momentum, energy, and entropy equations, the simplified one-dimensional 


relations 

are: 



d 

dt J 

f pvdV 

dU 

+ m ... = 

' v dt 

n k )rfi k - Y.P k n k A 

k k 

— 1 

dtJ 

f [p c p T 0 

V 

- P)dV = ll] 

L ^ ^ k{hu ' n A ) * // k r * 1 k ( 1 

k k 

-J 

dtJ 

f psdV 

V 

- Y.srn{i u -n) 

k 

> 0 


An important goal in transient ejector analysis is finding an accurate approximation to 
the time dependent volume integrals on the left hand side of each equation. 

Remarks on Supplementary Equations 

Overview 

Application of the general, one-dimensional flow, control-volume equations to 
the ejector of figure 10 is intended to provide a framework for ejector performance 
prediction in terms of the primary nozzle flow, initial conditions from which the 
secondary fluid is drawn, and ambient conditions. In application to the mixing 
region, the assumptions involved are critical since it is within this region that the 
most complex physical phenomena occur. 

In anticipation of the section 6 discussion on methods of solution for the 
complete (time dependent) one-dimensional ejector equations developed above, the 
simplification of steady-flow is explored in section 5. That excursion is of value in 
situations where a quasi-steady flow can be assummed to exist. That discussion also 
addresses some concepts to closure common to the unsteady flow problem 
formulation. This allows section 6 to focus on the difficult issue of evaluating the 
time-dependent volume integral terms. A summary of the solution options that exist 
and the reccommended approach in the provision of a complete system of equations 
is then made. 

The control volumes used in the applications to follow are illustrated in figures 
1 and 10. Since it is assumed (for this phase of the work) that the primary nozzle 
flow characteristics are known, no control volume is drawn for it. An important 
assumption is that the nozzle and the shroud are adiabatic - as discussed earlier this 
is not (immediately) an assumption that the flow is isentropic; the latter is allowable 
only if viscous dissipation is absent from the flow (as you would have, for instance, in 
an ideal, irrotational flow). 


- £ 


s dm > 0 
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To set the stage for application of the control volume equations to the ejector, 
some remarks on the character of some unknowns and the working medium follow: 


The state principle allows that local thermodynamic states are expressible in 
terms of two thermodynamic variables; in the present work the static pressure 
and temperature have been chosen. 


2 . 

3 . 

4 . 


5 . 

6 . 
7 . 


8 . 


The assumption of an ideal gas provides the entropy and enthalpy are 
functions of temperature only. Density is given as a function or I and r 
through the ideal gas law. 

The local mass flowrate is defined in the one-dimensional case as a function of 
velocity, temperature, and static pressure. 

It is assumed that the force on the shroud can be constructed as the sum of 
empirical relations for duct flows, including such effects as losses due to 
expansion or contraction of cross-sections. The empirical constants will reduce 
the unknowns in the function for F to velocity, temperature, and static 
pressure. Geometric characteristics arc assumed known. 

In general, entropy is employed in arguments related to the admissability of 
solutions, not in the direct solution for specific velocity, temperature, and 
static pressures. 

Velocities at the inlet duct and primary nozzle discharge are uniform and 
one-dimensional, but distributions internal to the mixing region are express- 
able in some "appropriate" self-similar form. 


The default unknown velocity is the relativ e velocity, v, since the absolute 
velocity, u, and the frame of reference velocity, U, are related by v - u - U. for 
the present work it is assumed that U is known. Actually U is generally not, but 
if the ejector GDE’s are coupled with those for the aircraft then U can, in 
principle, be determined (in general, as a function of time). 


Possible existance of shock waves is recognized, but an account here would 
reauire they are normal and stationary with respect to the shroud (diffuser, 
inlet, and mixing region) frame-of-reference. In the present work subsonic 
flight is assumed for the operational envelope the ejector participates in. 


We suspect in the primary flight mode the ejector is involved in will not be 
supersonic flight. 


Supplemental Equations 

The framework of analysis is constructed through application of the con- 
trol-volume equations to the inlet region, mixing region, and diffuser. The assump- 
tions of the previous section introduce the following supplemental equations: 


m = pv-nl 
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p = P/RT 
v - u - U_ 

F = F(P,T,v) 
h = /i(T) 

s = s(7) 


The principle unknows at this point are the following (field) variables: velocity, 
v(x,t)\ static pressure, P(x,t)\ temperature, T(x,t). The steady-state solution for these 
unknowns is discussed in the next section; Section 6 describes the approach for the 
transient situation. 
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5. STEADY FLOW ANALYSIS 

Overview 

Although it is the intent of the present work to focus on unsteady ejector flows, 
the steady-flow solution methodology is of interest for several reasons: 

The steady-state solution is the starting solution for the unsteady analysis, 

Many of the subroutines generated for the steady-state version of the program 
are common to the transient program; a successful steady-state solution is a 
useful check on those routines, 

The steady-state ejector performance program can be a useful theoretical tool 
for interpreting data from the NASA Lewis PLF. 

Presentation of typical steady-state ejector performance equations highlights the 
mathematical benefits of disposing of time-dependent terms in the ejector analysis and, 
by default, provides the system of equations one would use in a quasi-steady flow 
analysis. For instance, the section 6 assumption of a quasi-steady inlet and diffuser flow 
means the steady-state inlet and diffuser equations discussed below are applicable. 
Additionally, many details associated with the execution of the general ejector 
"problem" (and coincident with an unsteady flow analysis) are easily illustrated with a 
aescription of the solution of the steady flow equations. The present work draws 

f iartially upon the discussions of Belivaqua [1974], Salter [1975], and Alperin and Wu 
1983a, b]. 

System of Equations 

In a steady-state ejector flow scenario, time-dependent control volume terms 
would not be involved in the system of equations. Recalling the general system and 
auxiliary equations presented in section 4, it is important to note the steady flow 
assumption leaves three governing equations and (only) nine unknowns for the 
description of the mixing region. Although there is some variation in specific analyses, 
closure to this problem is usually obtained by prescribing one or more of the following: 

1. confluence static pressure ratio at the inlet, 

2. primary and secondary inlet velocities, 

3. static pressure at exit, 

4. stagnation temperatures at the inlet. 

Often, if there is not a direct prescription of, for instance, the stagnation conditions at 
the inlet, an isentropic flow is assumed between the (assumed known) free-stream 
conditions and the inlet cross-section, thereby facilitating calculations of the inlet 
temperature and pressure. 
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For the present steady-flow analysis, the three conditions for closure are obtained 
by prescribing (a) primary flow inlet conditions, (b) static pressure at the diffuser exit, 
and (c) an isentropic inlet and diffuser. A detailed discussion of the system of equations 
follows. 


Mixing Region 

Considerable simplification of the system of equations is realized by invoking 
the steady-state assumption; the system is given by: 

Mass: 

0 = m , P + rh ls - rh 3 ( 1 ) 


Momentum: 

0 = v lP rfi ]/ > + Vi S rfi {S -p<v 3 >rh :} + (f :> ]P -P 3 )A 3 + F 


( 2 ) 


Energy: 

0 = rh. | P // o , i p + m 1 5 // o , 1 5 — rn 3 // o , 3 


( 3 ) 


The entropy equation remains a condition on the solution, as discussed earlier (also, 
see Appendix B). Some assumptions often applied to the system are: 

1. All skin friction and blockage losses are neglected. 

2. The primary flow is fully expanded at the point of confluence of the two jets 
(equal static pressure at the inlet region, station 1). 

3. Specific heat of each gas stream is assumed equal. 

4. The ejector is stationary in space (test stand set-up). 

5. Gases are assumed thermodynamically perfect (constant specific heat) and 
described by the equation of state. 

Only assumptions 1, 3, and 5 are used in the present work. 

Non-dimensionalization of the system of equations above clarifies the partici- 
pation of each equation in the solution; related manipulations are outlined below. 

i. Continuity 

A useful non-dimensional form of the mass flowrate equation is (see 
Appendix G.6) 
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m 

~PA 


( RT ) 


1/2 


y| 1 + M 2 


1 / 2 


/ 6 ( y.A/) 


( i) 


where it is understood that the specific heat ratio is known or can be computed 
from Dalton’s law (Appendix F). From conservation of mass 

rhp + rh s - rh 3 = 0 (5) 


in non-dimensional form 

m 3 m p m $ n i p 

= + = ( 1 4 ) 

P P A P PpAp PpAp P f Ap J 
= ( p p r o . p ) / 6 ( y p • p ) ( i + p ) 

If the mass flowrate equation for station 3 is normalized by the primary flow 
conditions 

^3 _ P 3 A 3 f J , M 3^ 

PpAp PpAp y[R^f' 0 ~ (7) 

then by comparison with the result derived above 

P 3 _ j~^-3T 0,3 f 6{y P < ^ P ) 

Vp = T 3 y /?,7- 0 .,/ 6 (y 3 ,M 3 j { /i> ( } 

the entrainment ratio itself can be expressed as a function of in a similar form: 

^ _ rhs _ P s 4 s j~P p P o ./>/*( y s ’ M s ) 

rhp P p A /> y P sT o .s f ${y p , M p) 

so that the static pressure ratio is given hy 
p p _ a s j~p^p i i) i 

^ V PsT o.s f (>{y p,M p)p 

In general, the stagnation pressures are assumed to be related isentropically to 
the static pressures at the same station by the function 
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P 0 = i 1 + P = f ^ y ’ M)P 

then the stagnation pressure ratio is 

P 0,S _ P S I 3 (y S > M Jf) 

P o , p P pfAyp’M r) 

ii. Momentum Equation 

In non-dimensional form, the momentum equation is given by 

P S A S P 3^3 P P V P P pt) s P 3 < u 3 > 

1 + + + R = 0 

P p A p P p A p P p P p P p 

The velocity terms can be written 


so that the momentum equation can be written 


{l +y P M 2 F ) 


P s A s 
P p A p 


0 + y sMl) 


P 3 A 3 

P p A p 


{\+Py*Ml) = 


Re-arrange the previous result for the static pressure ratio 

P 5^5 _ lR S T 5.0 / 6 ( y P ' M p) 

P p A P V P p I p .0 / e(y S • M y] 


so the momentum equation is 


0 = (l + y P M 2 p) 


RsT 

R P T 


5,0 


P, 0 


l / 2 


/ 6 (y p.M f ) 
f Ays, Ms) 


/ 4 1 + y s M s) 


- 0-/0 


P 3 ] 3.0 
p P I P. 0 


/ 2 


ft{y_p,Mp) 

feiy-j'Ma) 


(\+py 3 Ml) 


then 


( 11 ) 


( 12 ) 


( 13 ) 


( 14 ) 


( 16 ) 


( 17 ) 
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0 


= /,( V,.M l 


R*T 


S 1 S 0 


R pT p 0 


i / 2 


f ?{y s * ^ s) 


~ ( 1+/0 


ft-,7 


I / 2 


3 7 3,0 


K#>T 


/> 1 p,o . 




(18) 


where 


/ 7 (y,A/) = (l +/?yAf 2 )// 6 (y.A/) 


(19) 


The solution for M 3 is now given in the form 


1/2 


/ 7 (y 3 .A/ 3 ) = ( 1 + /0‘‘ {fy{y P ,M P ) 

l R P * P, 0 


ft 


R S I s , 0 
R p I p o 




( 20 ) 


For given primary and secondary conditions the right-hand-side is a constant 

/ 7 (y 3 ,Af 3 )= K x (21) 

An expression for M 3 is derived by re-introducing the definition of the function 

Si 


1 + py M 2 = 

A' , M^y| 1 

(22) 

since 



( 1 + py M 2 ) 2 = 

1 + 2py M 2 + p 2 y 2 M 4 

(23) 


then 

1 + 2pyM 2 + p 2 y 2 M 4 = K] M : y + A ' 2 M 4 y( j (24) 

This is re-arranged to yield 
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1 +{2py-K]y)M 2 + [p>y 7 -K]y[? 7 


M' = 0 


SO 


/I - />V-/r?y|*~ 


(2i>) 


(26) 


S = 2{3y-K]y 


C = 1 


(27) 


(28) 


The solution for A /3 is therefore obtained from the solution to the quadratic 
equation 


a{mI) 2 + BM l + c 

which is given by 

-B±7fi 2 -4/)C 


0 


M: 


2A 


The root is expressed in functional form as 

^3 = f 17 ( /? . y 3 . K | ) 

///. Energy Equation 

Introduction of the entrainment ratio in the energy equation yields: 


m f _ m s _ 

^3^3,0 — CpT po + — CjTj 0 

m 3 m 3 


where c is the specific heat at constant pressure. Rearranging, 


m , 


3,0 


7 ' P 0 + /z — 7" 0 
mr + m. s [c 3 ’ c 3 *' 


(29) 


(30) 


(31) 


(32) 


(33) 


and therefore 
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1 cp 


( 34 ) 


T 


3,0 


( 1 + /')C 3 


r F , 0 + 



Expressed as a ratio: 

T 3 ,o _ 1 Cp I j + CjT j 

T p k o 1 + //c 3 \ c p7 p ,, J 


( 35 ) 


/V. Inclusion of Frictional Effects 

Salter [1975] provides corrections for friction effects through an equivalent 
pressure loss: 

AP, = (fLpv 7 ) / (2 D h ) ( 36 ) 


where / is the friction factor and D/, is the hydraulic diameter. Also, a 
momentum correction factor, K m , is applied to the mixing region exit to account 
for incomplete mixing of the primary and secondary streams; use of the 
momentum correction factor means the momentum equation is to be written 
explicitly in v (rather than in terms of the mass flow rate), the result: 

0 = p A ^ F u x F + p A | s v | s - K m p A 3 v 3 

+ (/>, ,-/»:, )>l 3 - AP,A 3 (37) 

Bevilaqua [1974] and Kentfield [1978J take a slightly different tack and invoke an 
incompressible flow assumption to simplify their analysis (they also remark on 
the error introduced in this approximation). 

v. Mixing Region Equation Summary 

Solutions for the pressure, velocity, and temperature at station 3 are 
obtained by analysis of the mixing region; since conditions at the primary station 
are given, the field variables at station 3 are functions of the unknown conditions 
at station IS. 

Manipulation of the equation of continuity in conjunction with the 
equation oi state yields the ratio of pressure at the mixing region exit and at the 
primary nozzle discharge: 

F 3 _ A P I /? 3 T o,3 f(,{y P ,M |>) , 

F p 4 3 y R P T o F / 6 ( y 3 , M 3 ) 

this is implicitly a function of conditions at station IS, via the entrainment ratio, 
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/«( yp’M P ) 

From the momentum equation, the Mach number at station 3 was derived in 
terms of the Mach number at the primary and secondary flow conditions at 
station 1, 

y 3 . k 3 ) 

Through the energy equation the relation between the stagnation temperature at 
station 3 and the stagnation temperature at the primary and secondary stations 
was found to be: 

P3.0 1 c p I £ s I j,o 

< 1 + n 

T p.o ( 1 + /i)c 3 V CpI p . 0 

In summary, the prediction of P3, V3, and 73 is a function of Ps, ^SO > 
and the pressure matching condition at the diffuser exit (more discussion on this 
point later). Analysis of the inlet region will exchange the unknown pressure, 
temperature, and velocity at station 1 for those at infinity. 

Recall that, for this particular predictive situation, the secondary inlet 
conditions are not known in advance since the secondary flow conditions are a 
consequence of conditions in the mixing region, not a cause for the same. No real 
reduction in the number of unknows is gained from the addition of these 
relations but is a necessary development for the complete system of equations. 
Relations for the diffuser are similar to those for the inlet and provide closure 
for the final system of equations. 

Here, it is of benefit to realize that for the thrust augmenting ejector the 
static pressure at station 4 should, in general, be equal to that of the ambient. 

The system of equations for the inlet and diffuser are discussed next. 

Inlet 

The inlet is given as the region between infinity and station IS; application of 
the principles of mass, momentum, and energy to that region are described below. 
As a first approximation isentropic conditions are assumed to prevail. 

L Mass 

For a single-stream flow, the mass flowrate is unity. 


rh p 
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so that 


P S V S A S = PnVmA 1 - 

or (see Figure 1) 

p s v s A s — p 0 v o A o 


(39) 


(40) 


since the characteristic area at infinity is undefined. 
ii. Momentum 

A balance of momentum for an ideal nozzle flow provides 

P S A s + m s v s = P „A„ + rh„v m (41) 


then 


p s / rn s v s \ _ A^(^m z i„ \ 

P~\ P s A s J ^J'V P <*, A „J 


(42) 


P S A S _ ( I 

P-A 1 - \ 1 +y s Ms j 


(43) 


An alternate expression for this ratio is given in terms of the mass flow function, 

/6 


PsAs m I RsTq.s \ 1/2 /6{ y -.-' tf -) 

P -Ac \R-To.-i fbiys’Ms) ^ 

but for an isentropic flow it can be shown the two ratios are not independent. 
iii. Energy 

Conservation of energy provides the simple statement that the stagnation 
temperatures in the free-stream and at station IS are equal 

T o.- = T 0 s (45) 

this result relies on the adiabatic inlet assumption. If local properties are related 
to local stagnation conditions isentropically, then 
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This is a slightly more convenient expression than that derived from the 
momentum equation since the area ratio is not required in the present result 
(but an isentropic flow assumption is). 

iv. Remarks 

In the general case of analysis the free-stream velocity, static pressure, and 
static temperature are known so the stagnation temperature and pressure can 
therefore be computed. If isentropic flow conditions are assumed to exist 
(adiabatic ejector, frictionless flow) then the stagnation temperature and 
pressure at station IS are equal to those in the free-stream. When the velocity at 
station IS has been determined, the static temperature and pressure can be 
computed. 

For situations where the isentropic flow assumption is not valid, we return 
to the momentum equation and account for frictional effects. As a first 
approximation of the "answer", an incompressible flow is assumed to exist. 
Recognize that in such a case, the differential formulation of momentum is 
preferred since Bernoullis’ equation can be extracted. This provides an equation 
applicable to any inlet streamline. The analysis of Salter[ 1975] introduces a 
correction factor for real-flow effects with the use of an inlet loss coefficient, K\, 
and the loss due to the impedance of the nozzles to the free-stream flow, Cq, 


P 


I s 




K 
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Given the free-stream velocity, conservation of mass is used to find the first 
approximation to the secondary inlet velocity. 



and therefore 
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Diffuser 

Pressure recovery in the diffuser is a function of the area ratio, and becomes an 
important parameter in ejector design. For the thrust augmenting ejector the exit 
static pressure is (generally) intended to he that of the ambient, but real flow effects 
may dictate otherwise. Nonetheless, in principle the overall (isolated) ejector 
pressure ratio is unity. 

Equations governing the conservation of mass, momentum, and energy in the 
diffuser are similar to those for the inlet, but a brief discussion is presented below to 
highlight the pressure matching condition at the diffuser discharge. 
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i. Mass 

For a single-stream flow, the mass flow rate is unity, 

^3 j 

rh 4 


( 57 ) 


ii. Momentum 

In the absence of frictional effects, the adiabatic diffuser is characterised by 


P<A a = / 1 + a \ 

P 3 A 3 \ 1 + y 4 A/ 4 / 


(58) 


In this result the Mach number is modified by the mixing effectiveness parameter 
when an account of non-uniform mixing is desired. 

Hi. Energy 

As before, the adiabatic diffuser assumption allows that 

T _A = (59) 

t 3 / 2 (y4-A/ 4 ) 


from which the additional assumption of isentropic flow permits 

P 4 _ /a (y 3 ■ M_ 3 ) 

p 3 / 3 (y4.A^ 4 ) 


(60) 


and since 


P o,4 P o . : 


(61) 


P 4 =P„ 


then the exit velocity can be computed from 


/ 3 ( y 4 . a/ 4 ) 



(62) 
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iv. Remarks 

For the assumption that the discharge pressure is equal to that of the ambient we 
find that the momentum equation can be used to predict the "appropriate" exit 
velocity for the given station 3 conditions. Recognize here that the station 3 
condition reflects on the original velocity assumption at station IS. The correct 
assumption will be that for which the mass flow rates at stations 3 and 4 are 
equal. 

As a first approximation to the compressibje flow solution the use of an 
incompressible assumption is of convenience; Bevilaqua [1975] and Salter [1975] 
providfe more detials on those approximations. 

Primary nozzle conditions 

All of the thermodynamic and flow variables required to describe the primary 
nozzle discharge can be derived from three input conditions, mass flowrate, 
stagnation temperature, and static pressure. A discussion of the relationships 
involved - including an account for choked flow - is given below. 

i. Primary nozzle pressure ratio 

Manipulation of the mass continuity equation provides the relation 


rh 

~A 



( 63 ) 


where 


P _ 0 
P 



M 2 


y/y- 1 


( 64 ) 


It is intended the nozzle pressure ratio be expressed in terms of the mass 
flowrate, stagnation temperature, specific heat ratio, and nozzle geometry. To do 
this, the equation above is rearranged so that 



Since the isentropic pressure ratio gives the Mach number as 


M 2 


2 

(y-i) 
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then 
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2(y- l/y) 


(67) 
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The nozzle pressure ratio is extracted from the solution to the quadratic 
equation; the result is: 
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( 68 ) 


from which we establish the relation (notice the "+") 


2 \ V/ 32 /l 2 yV 2 


i /2 \ y / y- 1 
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(69) 


This equation is valid for choked and unchoked flow; in the latter situation, 
however, there arises a convenient relation from the concept that the mass 
flowrate per unit area is a maximum for choked flow and coincides with a Mach 
number of unity. 

A solution for the nozzle static pressure as a function of stagnation 
pressure, stagnation temperature, and specific heat is given in Appendix J. 

Returning to the first equation of this section and eliminating the static 
pressure from the RHS provides 
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then for a Mach number of unity 
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For air at standard conditions 


y = 1 .4 


ft - lb, 

R = 53.304 - 

lb m R 


we substitute and obtain Fliegners’s formula 
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= 0.532 


P o 



The nozzle pressure ratio is then 




1 .88 
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which is in a convenient form for hand calculations. 
ii. Primary nozzle thrust for isentropic flow 

From the steady-state momentum equation the thrust of a nozzle, choked 
or unchoked flow, provides 

T = m t v e + A e { P o~ p <°) ( 74 ) 


Introduce now the idea for maximum thrust in which the area ratio is such that 
the exit pressure is the same as that of the ambient; furthermore, allow the exit 
velocity to be given by the relationship for isentropic flow. 


v. = 



( 75 ) 


The primary isentropic thrust is therefore given as 


T, = 


m e JyRT 0 


2 

(y- 1) 



( 76 ) 


It seems that for unchoked flow the static exit pressure is equal to the ambient 
(for the ideal case); if the nozzle pressure ratio is defined by 


N PR 



( 77 ) 


then the isentropic thrust equation takes on a convenient form for evaluation of 
flow conditions. 

For choked flow from a convergent nozzle the exit pressure is not, in 
general, equal to the ambient. It is then necessary in the analysis of choked flow 
to return to the original expression for thrust and account for this fact: 


< 37 > 



T s = rh a . yRT 0 


(y- 1) 


+ A a {P'-r a , 


( 78 ) 


where the subscript a referes to the ambient condition. Now, with the substitu- 
tions 

M. = P,v t A„ 

P . = P.RT. 

rn t rT o 


and the additional relations that 


- \j 2 c P I o \ I 1 


then the isentropic thrust equation becomes 
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Now, 
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T. lRT 0 {y- 1 ) 

T oV y 
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and for choked flow 
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then 
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so 
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Continuing, 

T. 
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which provides the result 

T = m/P^ 

Combine the terms 
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then we obtain the desired result for the isentropic thrust for a choked nozzle: 


( 88 ) 
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T = liiJTo 


(89) 


2R 

y(y + i ) 

Thrust, thrust augmentation, and ejector ejficiency 

As previously mentioned the ratio of the total momentum increment and the 
thrust that would be obtained from the primary nozzle under ideal conditions yields 
the thrust augmentation ratio, 

T SYS 

<t> = (90) 

r | P .IDEAL 


1 - 




where the system thrust is given by 

T sys - di 4 ( v A — i»„) 


(91 ) 


and the thrust from the nozzle under isentropic flow conditions is 

T | P.IDEAL = ^ I P V \P 


(92) 


where 
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The efficiency of the energy transfer process is given by the ratio of the kinetic 
energy of the efflux and the energy input at the primary nozzle; the energy efflux is 


KE . 


ni 4 /? < v A > 


(94) 


and the energy input is 

Tip = v t p A t p{ T o . \ p ~ T \ p) (95) 

Solution Options 

In the previous section there were (for the general case of analysis) presented 9 
equations describing ejector physics in terms of 12 unknowns. The difference between 
control-volume analyses generally centers around the manner in which closure is 
provided. Two basic approaches are distinguished in the present work according to 
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intent - the "direct solution" is used for parametric analyses and the "iterative solution" 
used for predictive analyses. The iterative approach is obviously preferred in the 
present work. 

Direct Solution 

From a computational point of view, a parametric analysis involves a direct 
solution to the system of equations. That is, it is a simple matter to specify a broad 
range of secondary inlet conditions and compute a corresponding set of discharge 
conditions. In those cases where multi-valued solutions exist, selection of the 
"correct" solution is assisted by the entropy condition. No attention need be paid to 
extracting which of the solutions are naturally occuring, only that the solution 
correspond with the prescribed conditions at the secondary inlet that, for whatever 
reason, in fact have arisen. 

The parametric approach is not in line with the objectives of the present work 
so there is no need for further discussion. Complete coverage of the parametric 
viewpoint is given in the works of Addy, Dutton, Mikkelson, and co-workers (see 
references). 

Iterative Solution 

Although the iterative solution involves the same set of equations as the direct 
solution, it is the objective in the former to seek the appropriate inlet velocity which 
provides an exit static pressure corresponding with the ambient. So rather than 
select arbitrary inlet velocities, a reasonable estimate at the correct value is used 
(through, for instance an incompressible flow analysis) that meets ambient pressure 
conditions, and the result then refined - by iteration - until convergence is reached. 
A typical solution approach is presented below, although, as alluded to earlier, other 
solution procedures given in the literature may vary in some details of the steps 
involved. 


i. A Typical Solution Procedure 

The procedure outlined below reflects the assumption of insentropic inlet 
and diffuser physics; the impedance of the nozzle to the flow and other 
"real-flow" effects are ignored for the present discussion. We have: 

1. Compute the free-stream stagnation pressure and temperature Use the given 
free-stream static pressure, static temperature and velocity in the following: 

V „ = | U_ | sin a 

M. = VJc„ 

P 0 . m = / > ./ 3 (y..A/.) 
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T 


0 , W 


2. Compute the primary nozzle discharge conditions Use the (prescribed) 
discharge mass flowrate, stagnation temperature, and static pressure in the 
following procedure (Appendix J describes the modification if the stagna- 
tion pressure is given instead): 

(i) Compute the NPR based on the assumption that the discharge pressure 
is the same as the ambient: 


NPR 



1 + 4 


rn 2 T o R ( y- 1 
P 2 A 2 y l~T~ 



1 


(ii) Compute the NPR for choked flow (M= 1) 


NPR' 



_Y 

y-i 


(iii) If npr > NPR' , the nozzle flow is choked; use Fliegner’s formula to 
compute the stagnation discharge pressure, 



and the actual NPR is then based on the ratio of this stagnation pressure 
and the actual back pressure. 

(iv) If npr < npr • , the nozzle flow is subsonic. Use the static pressure at the 
nozzle discharge (an input condition) to compute the ratio 


= fl + l -l l+ J ln R ( IZl 

Pxp 1^2 2 \ [p' 2 A 2 y{ 2 


1 /2 


y/y- i 


then compute the actual NPR from the known back pressure condition. 

(v) With the correct stagnation pressure now established, the computation 
of the static temperature and exit velocity can be made 

~P \ p ~ T 0 iir / f 2 {y Xf ,M i ? ) 

V \ p = M j fC | P 
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3. Assume MIS The solution to the steady-state problem is an iterative 
procedure; in the present work a "starting solution" must be given. An 
initial Mach number assumption of 0.01 has been found a convenient 
starting point for the computations of the present work. 

4. Compute station IS conditions. The known free-stream flow conditions, the 
Mach number at station IS, and the assumption of an isentropic inlet are 
sufficient to estimate station IS conditions. The isentropic assumption 
allows one to equate the free-stream and station IS stagnation temperature 
and pressure: 


T 0. 1 s 

P 0 ,« 

^ 0,15 

* To,- 

Then 


T 1 s = 

t o,\s^ / 3 ( y \ s * m 15) 

T ls = 

T 0 .ij// 2 (y 1 s' ^ 15) 


The mass flowrate at station IS tan then be computed; enough information 
now exists to compute the mass e ntrainment ratio. 

5. Compute properties at station 3. Apply the system of equations from section 
4.2. 1.5: 

T 3 , o = F 8 ( d ' T s .0 ’ 3 p .0 ) 

m 3 = / 17(^3 ’ y 3 ’ f 3) 

T 3 = T 0 3 / / 2 ( y 3 , M 3 ) 

P 3 _ A p [~R 3 7~ 0.3 f 6{y p < M p) | + 

P p A 3 \ R p T 0 , p / 6 ( y 3 , M 3 j 

P 3.0 = P 3 f siy 3 ’ M 3) 

The mass flowrate at station 3 can be computed with these results. 
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6. Compute properties at station 4 The isentropic flow approximation and the 
assumption that the discharge pressure is equal to that of the ambient 
introduces 

P 4 = P m 



This velocity represents that which, if the station IS velocity were correct, 
would match the stagnation condition at station 3 to the static condition at 
station 4. Once the static temperature is computed 

P 4 = P 0 , 4 f f 2 ( y 4 * ^ 4 ) 

enough information exists to compute the mass flowrate at station 4. 

7. Compare the mass flowrates at stations 3 and 4. In general it is unlikely the 
initial guess for the station IS velocity is correct, so the mass flow rate at 
station 3 will not be equal to the mass flowrate at station 4 (the latter based 
on a static discharge pressure equal to the ambient). If this is the case, the 

"surge" is given by, 

, * , , 
m “ m ^ m ^ 

and the mass flowrate at station IS given by 
1 . . 

m I s ~ 2" 1 

from which the new inlet velocity at station IS can then be computed. 

8. Repeat steps 4 - 7 until the solution converges. The convergence criteria is 
that 

| rh * | < e = 0.1 

9. Proceed with the computation of thrust, thrust augmentation ratio, and ejector 
efficiency. 

This solution procedure is summarized in figure 11. 

REMARKS: 
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1. Step 6 is based on an unknown value of the exit pressure (station 4), 
allowing for the general situation where the exit pressure is not equal to 
the ambient. Some modification of the procedure is allowable when the 
latter is true. In such a case the computed P 4 is compared with the ambient 
and modification of the inlet velocity based on the error in pressure 
prediction. 

2. There is required a second iteration contained within each timestep if an 
ideal gas is used instead of a perfect gas. Figure 12 provides the 
computational steps for the inlet region. These steps are typical for each 
control volume and differ only by the subscripts of the variables used. 

Sample Computations 

A FORTRAN program for the Steady-state Ejector Analysis (SEA) performs 
computations for the methodology outlined above (Preparation of a user’s manual for 
this routine is currently in-progress). 

For the purpose of comparing theoretical predictions with experimental data, 
some preliminary data from the DeHaviland ejector tests at the NASA Lewis PLF 
facility have been made available to the present work; an overview of this data is 
presented in Appendix K. 

Results 

Input data to simulate PLF ejector runs 223-239 were considered adequate for 
the test of the steady-flow methodology since 

A. The nozzle mass flowrates were between 18.7 lbm/s and 46.97 lbm/s -- a 

range broad enough for the primary nozzle to be choked somewhere 
inbetween (about 30 lbm/s), and 

B. The primary nozzle stagnation temperature of 760 °R was considerably 

higher than the 540 °R ambient stagnation temperature. 

Figure 13 illustrates that the steady flow analysis presented here is capable of 
predicting both the magnitude and the trends obtained in the experiments; note the 
mirroring of the "dip" in the vicinity of choked primary nozzle flow. 
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6. UNSTEADY FLOW ANALYSIS 


Focus on Transient Ejector Flow 

Fluid flows characterized by time-dependent velocity fields are termed unsteady. 
More revealing descriptions are linked to time-asymptotic flow behavior. A transient 
flow is a ’temporary’ unsteady flow, associated with, for example, a change in ejector 
operation from one steady-state condition to another. Contrast this with oscillatory 
flows in which a periodic time-asymptotic flow character is exhibited. Ejectors utilizing 
pulsed primary nozzle flows are of the latter type. In the present simulation the focus on 
transient, not oscillatory, ejector phenomena descends from flight-critical aircraft flight 
control scenarios; an example would be transition to forward flight from vertical 
take-off. 

Effects the transient ejector simulation should capture are: 

1. Ejector response to the primary jet actuation, 

2. Momentary depletion of net thrust due to reallocation of engine fan air to the 
ejector, 

3. Feeder line delays (related to actuator transients). 

Upon integration of the ejector, engine, and airframe simulations, incipient vertical 
deceleration effects (due to thrust lag) can be quantified. 

Remarks on the Energy Exchange Process 

In connection with recent research on the coherent structure of turbulent flows, it 
appears that even steady-state ejector operation relies on (local) unsteady flow 
processes*. That is, in a coherent flow there are continuously deforming boundries 
whose motion can fin principle) be tracked - not time averaged as in classical theory - 
so that there is a distinctly traceable mechanism for work by pressure to be done and 
therefore energy (attributable to pressure) between streams to be exchanged. Addition- 
al energy exchange is provided by viscous shear forces, but these do not rely on the 
motion of the boundary. Since the pressure-exchange process (non-dissipative) is very 
short in duration relative to the (dissipative) process of mixing (viscous-dominated), 
there is potential that a pressure-exchange dominated process is likely to be more 
efficient than a viscous dominated process. This notion is supported by the marked 
increase in performance of (correctly designed) ejectors with pulsed primary flows in 
comparison with their steady-state counterparts. 

Since the literature reveals that most unsteady ejector studies concern pulsed laser 
operation, it is of interest to employ the conclusions from these studies in the present 
work (especially for work containing comparisons with experimental data). Such may 


1 This "unsteadiness" in a turbulent flow has, of course, been known for a long time, it has just traditionally 
been a practice to use time averaged quantities in, for instance, the Navicr-Stokcs equations. See Liepmann 
[1979], and Hussain [1981] for discussions of "classic" vs coherent turbulent flow modelling. 
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not be the case, however, since there is a distinct domination of pressure-exchange 
mechanisms in pulsed jets that cannot be assumed to remain in a steady-flow turbulent 
jet entrainment scenario for thrust augmenting ejectors^. 

The present work employs the concept of kinetic energy exchange between the 
primary and secondary streams for the transient secondary velocity prediction. There is 
a limitation to this approach in that an emprical correction is needed for closure. 
However, there is significant technical antecedant for the correction factor and its use 
follows from the work of Korst and Chow [1966]. It is important to remark that the 
creation of a kinetic energy balance for this purpose requires deletion of the 
corresponding kinetic energy balance terms from the energy equation for the control 
volume (what remains is the thermal energy equation). 

Three Levels of Approximation 

From a technical standpoint difficulties in real-time ejector simulation arise from 
evasion and compromise. By "evasion" it is meant we are in search of the answer 
without mandating recourse to solving the full unsteady Navier-Stokes equations; in 
terms of "compromise", there is a need to balance the level of approximation used (for 
which empirical calibration of the theory is then intensified) and the expectations for 
the accuracy of the simulation. Three levels of approximation can be considered: 

1. Quasi-Steady Flow: Assume that the characteristic time for changes in the forcing 
function (boundary conditions) is greater than the characteristic time for response 
of the fluid - that the fluid is "very agile" and therefore permits the steady-state 
equations of motion to be used at each instant in time. 

2. Characteristic Volume Approach: Allow that the mixing region can be partitioned 
into three characteristic volumes; one domain characterizing secondary flow 
effects, another for the primary nozzle flow, and one to characterize the mixed 
flow domain (the size of the control volumes are time dependent). 

3. Finite Volume Approach: Identify a finite number of control volumes of fixed 
size, partitioned only in the streamwise direction (the size of the control volumes 
are time independent). This approach introduces a relationship between the 
primary and secondary flows with the use of the self-similar profiles. 

Remarks 

In the analysis of a given unsteady flow problem, it is quite convenient 
mathematically if a quasi-steady formulation can be assumed valid. This leads to the 
use of steady flow equations in an unsteady flow analysis; at each instant in time the 
flow is assumed to instantly respond to boundary condition transients. A basic issue, 
however, is whether the characteristic time of the forcing function is the same order 
of magnitude as the relaxation time of the flow. 

There is no need to remark on the details of the system of equations for this 
approach since this is, by default, already given in the steady flow discussion. 


2 An interesting note here is that, at the other extreme, for pulsed flow there exists an entrainment even in 
a laminar flow tor those geometries where a secondary streaming (a viscous phenomena) is present. 
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In the case of external unsteady flows the Strouhal number (a characteristic 
parameter for the unsteady flow frequency) plays a significant role as the criterion 
by which to (or not to) invoke the use of a quasi-steady formulation (see Drummond 
[1985, 1986]). Although the Strouhal number is also a convenient characteristic 
unsteady flow parameter in internal flows, there appears no technical antecedents 
on which to base the prescription of a threshold for quasi-steady ejector fluid-dyn- 
amic operation. Only testing would allow this assumption to be made a-priori. 

An appropriate application of the characteristic volume approach is for flows 
that exhibit distinct flow regime characteristics, like, for instance, the inviscid 
interaction region discussed in section 3. Time lags for this system are generated 
through imposition of lag coefficients between field variables of the characteristic 
region; although this is a simple approach, it relies on accurate knowledge of the lag 
coefficients for the simulation to be accurate (read: analysis cannot be divorced 
from transient experimental data). This would be a new method in the approxima- 
tion of transient ejector performance, but is not explored further in the present 
work. 


The finite volume application of this work is also new in ejector analysis, but it 
is anticipated to be more accurate than the characteristic method since considerably 
more elements are employed; as expected, though, execution is likely to be more 
CPU intensive. Here, a hybrid approach is based on the work of Drummond [1985], 
and Seldner, Mihaloew, and Blaha [1972], Note that an important feature of the 
method is that steady-state data (not transient) can be used in the calibration of the 
method. It employs control volume elements of constant volume and assumes the 
mixing region flow is expressable in self-similar form. 

Finite Volume Approach 

It is evident that the quasi-steady and characteristic methods of analysis previously 
discussed really make no specific assumption about the nature of the flow in the mixing 
region. The finite volume approach attempts to overcome this problem. Indeed, a 
fundamental philosophical point to be made is that in exchange for the ability to more 
accurately depict conditions (than you would otherwise have) inside a given sub-region, 
a more complex form of the surface integrals must be accommodated. 

In the present work the inlet and diffuser regions are permitted to be represented 
by a quasi-steady approximation, and the finite volume method of analysis applied only 
to the mixing region. This is based on the perception of the inlet and diffuser physics to 
be driven by imposed pressure gradients, and the mixing region dominated by turbulent 
viscous interaction. Obviously, the notion is that the physics of the mixing region are the 
cause of a situation that yields the pressure gradient effect in the inlet and diffuser. 

The fundamental assumption in the finite volume approach is that the mixing 
region is divided into N sub-regions of known volume and whose individual velocity, 
temperature, and pressure can be given by characteristic quantities. A significant 
departure from previous discussions is that the characteristic quantities are not 
necessarily uniform over the sub-region, but do relate to characteristic distributions. 
These distributions relieve us from a specific treatment of, for instance, an inviscid 
interaction region, though by default such a phenomenon should be predicted within 
the domain of an accurate solution methodology for the problem. 
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Corresponding with the finite volume method of analysis is the need for specific 
statements about the characteristic distribution within each sub-region. A discussion of 
the profiles used in the present work follows an overview of the basic control volume 
equations for the mixing region. Then, application of the finite volume method to the 
continuity, momentum, and energy equations for the mixing regions is discussed. Lastly, 
the complete system of equations for the ejector are assembled and the proposed 
method of solution presented. 

Figure 14 illustrates the finite volume nomenclature and some typical elements. A 
virtual grid representation is given in figure 15. 

General system of equations 

The basic form of the control volume equations are: 


Mass 


df 

dt 


{ 


m 


) 



Momentum 


d 

dt 


X 


pvdV + m cv 


cW 

dt 




n ■ (p J)dA + F 


Energy 


-f 

dt J v 


PdV 



y Pv • ndA 


(y- 



PU ■ ndA 


( 1 ) 


( 2 ) 


( 3 ) 


Entropy 


d 

dt 


X 


psdV 


X 


s dm. > 0 


( 4 ) 


Again, the presence of the surface integrals (that in the steady flow analysis were 
simple sums of average quantities) reflects the idea that one cost of incorporating a 
simple turbulence model in the time dependent volume integrals is given by 
increased surface flux term complexity. 

Note the heat equation has replaced the general power equation for the energy 
equation (see derivations in Appendix A). Application of the self-similar profiles is 
in fact easier for the heat equation than it is for the power equation, but both are 
eventually needed in the analysis. As mentioned earlier, the kinetic energy 
components will be accounted for in a separate balance intended to predict the 
secondary flow magnitude; this will be detailed later. 
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Turbulent Jet Approximation 

A diagram of the present turbulent jet geometric characteristics is given in 
figure 16. Extending from the mixing region inlet plane there exists a potential-core 
region characterized by a fairly uniform centerline velocity, with no transverse 
component. This is distinguished from the mixed-flow region where the centerline 
velocity decay arises from momentum transport to the entrained fluid. 

This section extablishes the basic features of the potential- and mixed-flow 
regions for use in the finite-volume analysis. 

/. Jet spreading approximation 

Two jet angles are associated with the growth of the turbulent jet. One 
portrays the decay of the potential core region and the other bounds the outer jet 
boundary layer growth. Rectilinear profile assumptions for both have an analytic 
foundation and empirical verification. Although the outer jet expansion angle 
depends on the ejector pressure gradient and inlet velocity ratio, the inner jet 
expansion seems more exclusively a function of the latter. To eliminate the 
introduction of additional unknowns in the ejector problem formulation it is 
worthwhile to explore an analytic approximation for the inner boundary length, 
bj ; see the potential core region of figure 16. Below, some remarks on the outer 
jet expansion follow a discussion of the potential core approximation. 

ii. Potential core region length 

Characterization of the potential core region is given in the present work 
by "calibrating" an analytic model; calibration is done with ejector data. 
Abramovich [1963] derives the length of the initial region for 2-D planar 
co-flowing jets as: 

L _ t + r 

b 0 ~ c,(l-r)(0 .416 +” CL 1 347) ( 5 ) 

where the velocity entrainment ratio is 

r = v i s J v i p (6) 


and the empirical coefficient for free jets given by 

0.2 < c, < 0.3 


Alternate models are presented by Abramovich for multiple jets that link the 
increase in thickness of the jet proportionally to the intensity of turbulence in the 
stream (those models also assume a loss of single jet identity loss for multiple jet 
configurations). We find it convenient to use the expression above with a 
corrected constant 
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c 


0.4 


which is based on the experimental work of Bernal and Sarohia [1983]. 
Prediction of L/b() of 15.5 for a 2-D single jet corresponds fairly well with a value 
of 18 extracted from their plots of nondimensional centerline velocity. In this 
regard the work of Krothapalli et. al. [1980] is of interest since multiple 
rectangular free jets were considered and the dimensionless length estimated 
(from their plots) to be in the vicinity of 14-20. This is the basis for the earlier 
remark that, at least for the ejector configurations of interest, the rate of the 
potential core spreading is influencedmore so by inlet velocity than longitudinal 
pressure gradient. For the limited purpose of establishing the non-dimensional 
potential core jet length the free turbulent jet results are applicable. Bear in 
mind that unlike the free turbulent jet an axial pressure gradient is assumed for 
the ejector. 

Figure 17 illustrates the non-dimensional centerline velocity as a function 
of the non-dimensional centerline distance for a typical multiple free jet. The 
potential core region is important since the total non-dimensional length of the 
mixing region under consideration is on the order of 20-30. 

til Outer jet boundary expansion 

A linear representation for the outer jet boundary is illustrated by the 
boundary h// in figure 16. This is based on the interpretation of data by 
Abramovich for incompressible planar jets (free and submerged). Data from 
Donsi et al [1979] supports this trend for extremely large pressure gradients 
(experiments were for fluidized beds with jetting). We therefore have 


and note Abramovich’s simple expression for the constant has the attractive 
result 


b 

z 


c 


1 -r 
I + r 


( 8 ) 


where ci is the same as used previously. An account of jet growth as a function of 
longitudinal pressure gradient could be nested in the constant. Due to the 
time-dependence of the pressure gradient for transient ejector operation the 
(outer) jet growth -- and therefore the proportionality constant — is also a 
function of time. Since the development of Abramovich is based on continuity 
the present work follows suit, but from a finite volume standpoint. This ensures 
the finite-volume method will satisfy continuity without an ad-hoc construction of 
b(x). 


< 51 > 



It is interesting, though, to remark on a test case for the correlation above 
where r=0.35 and an ejector mixing region pressure gradient of 14 Ibf/fr is 
given. The 9° expansion angle from finite volume continuity considerations 
compares reasonably well with the 10.1° value from the Abramovich model 
(with C7=0.4) and, as expected, is higher than the 5° value found with cj for a 
free jet with no longitudinal pressure gradient. 

Self-similar profiles 

In the characterization of an element of the mixing region the velocity 
description simplifies with the 2-D planar turbulent jet self-similar profiles of 
Abramovitch [1963] for co-flowing jets: 



Alternate co-flowing jet profiles are used, for instance, in the work of Korst and 
Chow [1966] or Lund, Tavella, and Roberts [1966], but the polynomial form is more 
applicable to the physics of interest here. 

Introduce the dimensionless radial coordinate. 



so that in general form 


V - V, 
V m V q 


*(?) 


then 


/( 1 - e 1 5 ) 2 , 0 <e < 1 \ 
\ 0 , 1 <?<| / 


V =u # (l-0) + u m <t> 


/(*) 


( 11 ) 


( 12 ) 


where 

| = B/b 


(13) 


and B is the maximum value of the jet half-width, b. For an incremental volume of 
length Az, width y, and height x, only part of the jet consumes the element. 
Therefore, the general velocity profile is written 


v(f) 


//(*). 0<*<M 

\ v„ 1 <f <| / 


(14) 
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The description above makes no account for the potential core for the initial 
jet region; for this we note 


b = 


b n b i 


( 15 ) 


and the dimensionless radial coordinate becomes 

x - b , x - b , 

* b a ~ b , b 


( 16 ) 


Recognizing the inner boundary vanishes outside the potential core region (by 
definition), then by default the dimensionless coordinate takes the correct form if 
we define 


where z/ is the length of the potential core region. In the discussion to follow 
distinction between the potential core and mixing regions is not necessary if the 
appropriate non-dimensional parameters are understood. 

To allow for a density variation across the jet a basic self-similar profile is 
assumed, 


ztanO, z < 7. { 
0 ,z> z , 


( 17 ) 


P~ P. 

Pm P o 

where 

/l(H = 

and obviously 
P = P. 


/l(0 

fl, o<f<r\ 

\o, 

1-/1) + Pm A 


( 18 ) 


( 19 ) 


( 20 ) 


The value of t ' used in the present work is unity, though fluidized bed data suggests 
a value of 0.9 allows more interaction between the jet boundary layer and the 
free-stream. For a constant cross-sectional density there is no need for a self-similar 
approximation; thus 

p = p ( 21 ) 


Pressure is assumed to have a profile similar to that for density. 
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p- p 


A(Z) 


( 22 ) 


P m -P t 


Again, this profile need not be invoked if the assumption of a uniform transverse 
pressure gradient is specified. 

Temperature profiles are shown by Abramovich [1963] to have the dimension- 
less profile 


T - T 0 

T n ~T „ 


*(?) 


1 


i s 


( 23 ) 


derived as the square root of the velocity profile. 

Transformation of coordinates require that derivatives of f be determined. If it 
is assumed that b is in general only a function ofz and t (not*), then. 



so, for example, the generic area integral across the mixing region is, 


£ = 



z(p)W dx 


= W | j z(p)dx + J~ z(p)dxj 
= Wbl [ J\(p)dZ + J^z(p)df J 

and the integral for the complete cross-section would be 


Z 


I C.G’rt) 


i- 1 


( 25 ) 


( 26 ) 


/. Example Application: Average velocity computation 

Computation of the average mixing-region velocity at station A: is a simple 
example of the use of the self-similar profile integrals. Define the average 
velocity by 


v 


k 



vdA 


( 27 ) 
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we have for "AT nozzles that 


2NW r B , 
V k = — — vdx 

A k Jo 


and through a change of variable 
2NWb K ft 


A k 


v,dt 


Application of the self-similar profiles yields 


2 NWb k r t 


A, 


f {y.(l + 

J o 


where 


_ /( 1 - f ' s ) 2 .0<f < I \ 
1 \ O.I<?<? / 

Completion of the integral yields 


2 NWb k j n 

X,' (""Jo /l<<{ 


k. J ( I ( I )df 


Let 

2 NWb k 


then 

D = A*{u m ( 1 -0.8 + 0. 2b) «-w # (0.8-0.2 , o) + w.(f - I )} 


and finally 

v = yl ’ { 0 . 45 + - 0 . 48 ) } 

Application of the finite volume method 

The assumptions and profiles outlined in the previous section are applied 
equations of mass conservation, momentum conservation, and balance of energy. 


( 28 ) 


(29) 


(30) 


(31 ) 


(32) 


(33) 


(34) 


(38) 


to the 
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Application to conservation of mass 

For the generic sub-region k, bounded by surfaces at k and k + 1, 


dm A 
dt ) k 


= m , 


- rn 


*+ 1 


(36) 


where 


rn i 


2 N W I p k v k dx 


(37) 


From the similarity profiles 


rh k = 2NWb k J {p,(l - A) + p m A}{u,(l -4>) + i/ m <f>}df 


(38) 


As discussed earlier, an account of real-flow effects is provided in part through an 
empirically obtained jet spreading function b(z,t) that defines the jet boundary. This 
is a common function explored in turbulent jet analyses. In general, the spreading 
function is dependent on space and time, although for the present preliminary study 
considerable simplification is obtained by ignoring the time dependence. It suffices 
for the present work to admit a quasi-steady approximation for the jet expansion; 
this approximation can be checked when ejector data is available. 

After substitution of self-similar profiles the mass flux integral becomes 

rh = 2NWb + P e v m ( \ - A)<t> + p m v a A{\ - <t>) 

Jo 


+ Pd'„( 1 -/1)(I -*)W (39) 


then (see Appendix H) 


F ' " j 

f A<t>d% = 0.48 

(40a) 

■ J 

O 

II 

"tJ 

'N 

1 

(406) 

^ • J 

f/l(l-0)d^ = 0.88 

(40c) 


F< = J(l-/1)(l-*)df = l- 1 ( 40 d ) 
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where 


B 

b 


( 41 ) 


and the numerical values are for the mixed flow region. 

Now, 

m = 2N\Sb(p m v m F x + p t v m F ,+ p m v e F 3 + p e v t F 4 ) 

= 2NWb{Z x ) (42) 

The mass in the elemental volume is also given by 

m = f odV = 2NW Az f pdx 
J V ' ^0 

= 2NWbAz f pd{ (43) 


where the characteristic density approximation for the finite volume has been used 


P 



=4 P*. i 


(44) 


If now 



P n f l + P.^2 


where 

= f Adi- 

Jo 

E 2 = f\l-A)df 
J 0 


(45) 


(46a) 


( 46 b ) 


then 

m = ?.NWbAz(p m F. x +p.E 2 ) (4/) 

The continuity equation now has the form 
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(48) 


BAzl E 




dp e 


dt 



I k 1 Z 


) 


i. Uniform transverse density approximation 

A uniform jet density in the transverse direction simplifies computations by 
eliminating the distinction between entrained and primary flow density. This 
does not ignore the marked extreme of secondary amd primary flow, rather, 
treats the combined flow with a characteristic density extracted from thermal 
exchange upon mixing. With this assumption the elemental volume mass 
becomes 

m = 2NWAzpB (49) 


and the result for the finite volume density derivative is 

f dp\ _ b* Z i* ~ b*, , Z ( 

V dt J k BAz 


(50) 


ii. Computation of o„ 

Continuity conditions across the first mixing region finite volume (steady 
flow case) allow an approximation for the outer jet spreading angle. The 
approximation is based on the following 

1) constant secondary velocity, and 

2) constant primary and secondary flow densities. 

In practice, these assumptions provide the jet spreading angle as a simple 
transformation from "top-hat" profiles to those of a self-similar form. Since 


m , 

= 2NWb t Z , 

(51) 

then 



b , 

171 | 

2NW Z i 

(52) 

where 



b, 

= b/,-6, 

(53) 


and the outer boundary layer then defines the jet expansion angle. A more 
convenient expression is given by an explicit representation of b as a function of 
the primary and secondary inlet velocities. Recall that 


< 58 > 


Z, = p n v m F 1 + p'V m F 2 + p m o.F 3 + p'V'F 4 (54) 

where 

F , = — + 0.45 (55a) 

1 b 

F 2 = 0.0 (55b) 

f 3 = 0.55 (55c) 


By substitution 

m , = 2 Af 1/ ( p m w m ( b / + 0 . 45 b ) + p m u , ( 0 . 55 b ) + p , u , ( - b ; - b ) } (56) 


then the expression for the boundary layer thickness is 

m , - 2NV p{b ,v m + [B - b ,)v a ) 


b - 


2 N W p { 0 . 45 v m - 0 . 45 v „} 

The outer jet expansion angle is theref ore 

b , + b ; — b o 


(57) 


6,i = arclcin 


( b/, - b 0 ) 

_ [ 

J 

" l 


(58) 


The expansion angle is determined at the time of initialization of the flow and, as 
a first approximation, remains constant in time thereafter. This assumes the 
inner boundary layer is not a deciding influence on the overall jet boundary layer 
characteristic. The form of the inner boundary layer used in the present work is 
primarily a function of the primary and secondary velocity ratio. On the other 
hand, the outer jet boundary is influenced more by the longitudinal pressure 
gradient than by inlet velocity ratio. Providing that the free-stream (inlet) and 
discharge pressures remain fairly constant, the constant boundary layer assump- 
tion will, as a first approximation, not significantly impact the gross characteriza- 
tion of ejector physics. 

Application to the momentum equation 

The integral form of the momentum equation for the elemental volume is 
given by 
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(59) 


J t J v PUdV + m cv = -^pu(v-n)dA - j> n-(p/)dA + F 

i. Time rate of change of momentum 

Computation of momentum within an elemental volume is given by the 
volume integral of the density and velocity product, 

= J pvdV = 2NW Az J pvdx (60) 


For simplicity the average density concept from the continuity equation discus- 
sion is used here; furthermore, an average element velocity is assumed character- 
ized by the velocity at station k + 1, 


M 


k 


2NWAzp J v k cix 


(61) 


this yields 


M k = 2NWbAzp{v n [F y + F 2 ) + u q (f 3 + F a )} 


from which the time derivative is 

rf \4 . /-/ s-\ 

- 2W6Jz( 


dt 




(62) 


(63) 


ii. Momentum flux 

The momentum flux across station k is 


M , 


pv 2 dA 


= 2NW | p k v k dx 
J 0 


= 2NWb k p k v 2 k dt 


(64) 


Scalars are used in the above since the directions (signs) are understood at each 
station to be given by: 
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(65a) 

(656) 


“(**•«)* = +1 

Recall that 

v ■■ u.(l-*) + 

where 

(i-f 15 ) 2 . o<?<i\ 

0 , 1 <£<| / 

then 

u 2 = u?( 1 - <P) 2 + v 7 m <t > 7 + 2v a u m {<t> - <#> 2 ) 

and so 

M = 2NWb f (p,( 1 - /1) + p m /l){u 2 ( I - <t>) 2 + v 7 m <P 7 + 2.u a u m { 

J 0 

Performing the multiplications, 

M = 2NWb f <p.u 2 (l -0) 2 (1 -A) * p m u. 2 (l -0) 2 4 

%/o 

+ p.u 2 (l-/1)0 2 + /m u 2 /10 2 
+ 2p.u.u m (l -/I)!*-* 2 ) + 2p„w.« m /l(0- 

Introduce a set of integrals G/ 

C , = f A<P 2 d£ = 243/770 

C 2 = J/1( 1 -<*>) 2 df = 320/770 



-* 2 ))d* 

(67) 


2 )w 

( 68 ) 

(69a) 

(696) 
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c 3 = 2 j A{<t> - <P 2 )d% = 414/1540 

C 4 = f(l-A)<P 2 d{ = 0 
C 5 = J(1 -/l)d -0) 2 df = |-1 
c 6 = 2 J (l-/1)(<£-0 2 )d£ = 0 


(69c) 

(69d) 

(69c) 

(69/) 


so that 

A/ = 2/Vl/b{p m t/ 2 G, + p m v 2 C 7 + P m v e v m C 3 + p a v 2 n C 4 

+ P.^C S + P a V a V m ^ 6 } 

A/ = 2 /V 1/ 6 { Z 2 } 

The net momentum flux is therefore 

~f pv 2 {i v - n)i v dA = ( A/* - A/*, , )«„ 

- 2/VI*'(b,Z 2> -b,.,Z 2> „)(*i.) 

iii. Static pressure integral 

The surface pressure integral is re-written as 

-fP(l-n)d.A = - f P( + i • (-t))rM - f P( + i-i))dA 
J A J A * J I 

so application of the self-similar profiles is restricted to the equation 
f A p dA = 2NWb C Pdt 

from which we obtain 


(70) 

(71) 


(72) 


(73) 


(74) 
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L 


PdA 


2NWb j\p t {\-A)+P m A)dl 

2NWb{P m E x + P 0 E 2 ) 

2NWbZ 0 


where Ej and £2 have been previously defined. 

The surface pressure integral is now given by 

-£/>(/-n)cM = +2NW(b,Z 0 ' k -b k . i Z 0 ' k . l ) 


(75) 


(76) 


iv. Wall friction 

Frictional effects are considered an explicit function of velocity, pressure, 
and temperature (or density) 


F k 


Ek{p,’ v i>Pi) 


(77) 


but in the first approximation frictional effects are neglected. 
v. Summary of components 

The modified form of the momentum equation is given by summary of the 
components discussed above, the result is 

2NWb k . | J 2 {^ m (F 1 + F 2 > y ,(F 3+ F 4 )} 


+ -—p{E [ + E 2 ) + d>, r ° P {E 3 + E 4 )> 


dt ' v ' " (it 

= 2NW{b k Z 7 -b t , l Z 7 ) + 2Nls(b k Z 0 -b k . l Z 0k ) (78) 


The last term has the simple interpretation of the net momentum flux, 

M Nat = 2NWb k [z ik + Z 0t ) * 2 NWb k , x {z 2k ^Z 0t ) ( 79 ) 

where at the ejector inlet 

= ^ I P ( F 1 f + P ! P V | p ) + A , 5 ( P , s + p | S V | s ) (80 ) 
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A solution for the time derivative of velocity can now be extracted 


dv m _ A/ Wa( dpi ( + f / i 4 '] 

dt ~ 2 A/t/6** 1 ~Az dl p\ Vm + V \r , +F?)j dt\F\ + F~ 2 ) 


Application to the energy equation 

The form of the energy equation most conveniently adapted to the finite 
volume analysis is the heat equation, given in integral form as 


d 

dt 



-y 


pv-ndA-(y- 



■ ndA 


(82) 


Distinction between this interpretation of the energy equation and the mechanical 
and general power equations is described in appendicies A and B. Term by term 
evaluation follows; in all discussions use has been made of the uniform transverse 
pressure and density assumption. 

i. Surface flux terms 

Since pressure and density have similar representations, the pressure 
energy-flux term is similar to the mass flux term; by direct analogy 

e p = -yipv-ndA = d^-d*,, (83) 


where 

d* = y2NWb k p{v m {F l + F 2 ) + v a [F 3 + F 4 )} 

= y 2 /V 1/ fc> Z 4 (84) 


and where the self-similar profile integrals F/ are the same as defined previously. 
The pressure-energy flux associated with the free-stream is 

-(y-l )j pU-ndA = (y - \ )U A[p k - p k . , ) (83) 


If desired, the area term can be re-written as 

A k = 2 NWb k {F, + E 2 ) 


( 86 ) 



2NWb k $ k = 2NWB 


( 87 ) 


A k = 


ii Time rate-of -change of energy 

The volume integral of pressure cast in finite volume form is (see 
discussion on conservation of mass) 

J pdV - 2NWBAzp ( 88 ) 


and has the time derivative 

— f pdV - 2NWBAz ( ^ | 

dt J H \dt j k 


( 89 ) 


Considerable simplification of the time-derivative is realized from the approxi- 
mation that static properties at station k are characteristic of the finite-volume 
element. This is not a statement of uniform volume thermodynamic properties, 
rather, that the reference for derivative computations can be approximately 
given at either station k or k + 1. 

iii Summary of components 

In modified form the heat equation is 


2NWBAz 



y 2 NW(b k Z x -b k . x Z,) 

+ (y- i )t/2yvi/(b*z 0 -b^,z 0 ) 


from which 



-^(b k z x -b k . x y x ) + y -^u{b k z 0 -b k .,z Q ) 


( 90 ) 


( 91 ) 


Homentropic secondary flow 

Analysis of a finite-volume element whose cross-section spans the mixing 
region provides a relationship in time and space between the primary and secondary 
flows. The description of the mixing region is incomplete, however, since the basic 
profiles used previously only represent the characteristics of the jet boundary layer 
and not those of the secondary flow itself. In this situation an inviscid, homentropic 
flow assumption for the secondary flow provides closure for the mathematical 
representation of the mixing region. 
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The present analysis is predicated on a shock-free flow. However, as men- 
tioned by Anderson [1970a,b], introduction of shocks could be explicitly introduced 
through discontinuities in the initial flow field. 

Homentropic flow is distinguished by the absence of temporial and spatial 
gradients of entropy 

^ - Vs - 0 (92) 

Dt 


Entropy jumps across ei shock violate this condition. Nonetheless, in this work the 
fact that entropy is assumed uniform throughout allows the fluid flow to be 
described by the relation 

— ~ const (93) 

P y 

This replaces the need for the energy equation, but adherence to conservation of 
mass and momentum remains. The derivative of this expression yields 

dP_ _ Py 
dp p 


and thus 



(94) 


Conservation of mass allows for the computation of the density gradient in the 
secondary flow, but the representation chosen here is to compute the gradient one 
time step after changes in the primary flow have been computed. A point function 
computation arises 



At 


(95) 


A description of the entrained velocity derived either through a solution of the 
momentum equation for the secondary flow, or, as in the present work, use of a 
kinetic energy transfer function. The momentum equation approach is noted briefly 
below. Details on the kinetic energy function are given later. 
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L Momentum 

For an inviscid flow assumption Euler’s equation represents the momen- 
tum equation for a material volume. Although previous discussions consider 
extention to a material volume, the assumed existance of a streamline for the 
homentropic flow is a valid simplification. For one-dimensional flow we have 

dv dU , ...dv 1 dP 

— + — - -(u+t/)- — 

dt dt dz pdz 


Finite difference approximations for the streamwise derivatives give the time 
derivative of velocity in terms of velocity and pressure at stations k-1 , k, and k + 1. 

Finite volume initialization 

Nomenclature for the virtual grid used in the finite-volume initialization is given 
in figure 15. Initialization of pressure, velocity, and density is based on steady-state 
versions of the mass, momentum, and energy equations. The assumption is made that 
the uniform pressure and density profiles apply. Application is discussed below. 


Conservation of mass 

Elimination of time derivatives in the mass conservation equation produces 
0 - 2NWb t Z x -2NWb j Z ly 

where 

Z, = P m v m F l + p,v m F 7 + p m v,F 3 + P'V,F 4 


(96) 


(97) 


Now since 

m, - 2NWb , Z , j 
then the centerline density is 

P^m( /7 l +f 2) +1) .( /r 3 +/r <)) 

which becomes 

P,{w* + «.( F i + f . J} ' 2ATU by(F , + Fj) 

Solving for density, we obtain 


rn , 

2NWb, 


rrt, 


(98) 


(99) 


( 100 ) 
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a 


( 101 ) 


Pi 


v m + b 


where 


rh _ i 

2NWb j {F F 2 ) 


( 102 ) 



(103) 


Conservation of momentum 

Computation of the jet centerline velocity derives from the momentum 
equation and assumes the entraned velocity is known. First, recognize the discharge 
momentum is 


M | * A | />( P i p 


P i pu I F 


) + A 15>( p I ■ 


P \s v is, 


(104) 


and then at station k 

M lk = 2Wb t (z,/Zj 

where 

Z 0 “ {P m E^P.E 2 } 

Z 2 - { p m v 3 m G t + P m v]G 7 + P m v.v m G 3 

+ p ,v 2 m G A + p t v] G $ + p t v m G 6 } 


then 


A/, 

2NWb~j 


p{E 


£ 2 ) + p{u*(C l +C 4 )+u^(C 2 + C 5 ) + y.y m (C3-*-C 6 )} 


(105) 


(106) 


(107) 


(108) 


Normalization yields 
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M i 

2NWb i [G\ + G A ) 



f. 

G i 




Cj + Cj 

c, + g 4 



G 3 + G 
G i + G 


which has the shorthand notation 
p{v 2 m + v m c + d) + pe - J 


where 

c 


f Ga + G^ 

•Ig, + gJ 


d 


G 2 -G 5 A 

G , + G 4 J 


G i + G 4 
W , 

2Wb y (G, + g7) 



Conservation of energy 

The steady-state result for conservat ion of energy is 
0 - f pu ■ ncM 

= 2/Vl/{b*Z 4t -b ; Z 4/ } 

and therefore (by analogy with the conservation of mass discussion) 


P 


g 


+ b 


where 

*4. 

g 2/VU / b / (/ r i + F 2) 


)) 

(109) 

( 110 ) 

(111a) 

(111b) 

(111c) 

(Hid) 

( 112 ) 

(113) 

(114) 
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This transformation of the thermal energy exchange is supplemented by the constant 
enthalpy relation applicable to homentropic flows (see earlier discussion). 

Solution of the system of equations 

Substitute the density and pressure representations into the simplified momen- 
tum equation so that velocity remains the only unknown 


— — — [v 2 + v c + d} + — — — - / 
v m + b X m m ' Vm + b 1 


( 115 ) 


then by re-arrangement 


2 , <7G 

a 




CH6) 


/'Ld+S!-*/ 

a J a a 


0 


( 117 ) 


+ + C = 0 

This quadratic equation has the solution 

B 1 r-= 

*- ' -2‘7* - 4C 

where 

^Lif Lilli) 

till tc , + c 4 J 

^4,0 f If 1 + E 2 'j A/ | /^ Z 7 3 ^ 4 \ 

tii\ \C] + C 4 J V ‘rn\^rJ 


\C [ + G 4 


\G i + G 4 


( 118 ) 


( 119 ) 


( 120 ) 


( 121 ) 
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Kinetic Energy Exchange 

Analysis of the primary and secondary flow interaction has not, to this point, been 
completed. By themselves, the self-similar profiles close the loop for steady-state flows , 
but not transient ones ! This section provides an approximation for the turbulent flow 
kinetic energy exchange mechanism to characterize the influence of primary flow 
changes on tne secondary flow. 

Kinetic Energy Computation from Self- Similar Profiles 

Kinetic energy can be computed by the integral of the product of velocity and 
momentum (as described in Appendix B); from this there results the scalar quantity 

KE - - f pv 3 dS (122) 

2Js 

By substitution for the self-similar velocity profiles and with the use of a uniform 
density approximation, the kinetic energy becomes 

rt i 3 

KE = NWbp I [v e (] + d{ (123) 

•'f. 


If it is assumed 


f 2 > 1 


then 


(124) 


(125) 


KE - Wbp^%.( 1 -0) + tv^} 3 d? + 

Upon integration 

KE = NWbp{v]H x +v,v 7 n H 2 + b m v]ll 3 + v 3 m H < + u 3 ,H 5 ) 

where the Hi integrals are given in Appendix H. 

Change in Kinetic Energy of Secondary Stream due to Mixing 

Computation of the gain in secondary flow kinetic energy can be made by 
direct extension of the general kinetic energy relation above 


AKE - 


(126) 
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where ? defines the jet boundary streamline (for which the secondary mass flow 
through station i is equal to primary mass flow through /j. For the present discussion 
this dividing streamline is assumed known; Appendix L illustrates the typical 
approach of analysis. Expanding the equation for the change in kinetic energy yields 


AKE - Wbp|J^(u 3 -uu?)df + ^ (u 3 -uu 2 )df J 


(127) 


Consider the integration as the sum of the following four components; 




dt 


- v 2 H , + v t v 7 m H 2 + v^H 4 






KH 5 


( 128 a) 


(1285) 


X 


vv]dt 


v m v. F i +v , F 3 




d? 


3 r~* 

= 


(128 c) 


( 128d) 


In sum, the change in kinetic energy of the entrained flow is 

AKE = NU/6p< l ;„ 3 (W 1 + /y 5 - J F3-F 4 ) + u m u. 2 (//3-F 1 ) 

+ v 2 m v,H 2 + v 3 m H 4 } 


(129) 


Change in Primary Flow Kinetic Energy due to Mixing 

In a similar way as the change in secondary flow kinetic energy was computed, 
the energy loss of the primary flow is given by 


AKE 




dt 
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where the limi ts of integration reflect interest in the domain of the primary jet 
cross-section. 

Evaluation of the integral at station i yields the result 

AKE “ NWbp[v 2 H , + v t v 2 m H 2 + v m v]H 3 + v 2 m H 4 - v m F , - v,v m F 3 ) 


(131) 


Kinetic Energy Balance 

Computations for a specified steady-state condition show that the change in 
kinetic energy due to mixing is not the same for the secondary flow as it is for the 
Drimary Figure 18 illustrates some typical results (for net changes in kinetic energy 
between the* 1 input and discharge of the mixing region). In fact, the gain m kinetic 
energy of the secondary flow is entirely due to the mixing process, while the mixing 
loss of the primary flow is only a fraction of its total loss. In balance, however, the 
total change of kinetic energy of the primary flow is greater than that of the 

secondary flow. 

In the works of Korst and Chow [1966] and Chow and A ^dy [19M] the 
relationship between the change in entrained flow kinetic energy and the to 
primary flow kinetic energy is given by 


AKFjj_ a 

1 3 ^ 

2 P i p v i r 


(132) 


where a value of 12 for a for turbulent flow provides a reasonable match between 
theory and experiment for low-speed flows. At higher speeds the relationship 


12 + 2. 758 M 


1 P 


(133) 


is sometimes used. The important feature of this result is that the change in 
secondary flow kinetic energy has the functional form 


AKE , , 




(134) 


The difficulty with the energy transfer function described above is that it 
provides (by default) a quasi-steady flow approximation. It therefore cannot be used 
For the transient flow in its present form. To entertain local changes in pnmaiy-se- 
condary kinetic energy in a way that does not burden the computational procedure, 
consider the modified function 


AKE IS - P{AKE„. m ) 


(135) 
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where the subscript m denotes the loss in kinetic energy due to mixing. Numerical 
testing suggests that 


[AKE is ) 


t* M 


(zl/f 


+ C, 


o 2 z 


( 1 + M 


■UKE„ tm y 


(136) 


i r) 


This enhancement arises from the assumption that local velocity gradients in the 
steady flow case are typically less that the gradients experienced m the transient 
mode. Here, the introduction of an engineering approximation also results in the 
introduction of an undetermined constant, Cj. The alternative is to establish N 
computations of the kinetic energy exchange to coincide with the N control volumes 
of the mixing region; the present method permits post-processing KE information at 
the completion of mixing region calculations. Sample calculations indicate 

C, ~ 0.95-1.00 (137) 


Summary of Method 

This chapter has introduced the continuity, momentum, and energy equations 
required for the analysis of the ejector mixing region. A quasi-steady flow assumption 
for the inlet and diffuser relieves us from repeating the chapter 5 discussions for those 
components here. The energy equation is supplemented by the kinetic energy exchange 
equations presented above. Figure 19 summarises the solution methodology for the 
unsteady flow problem. As expected, several steps are common with the steady flow 
solution, especially in the initialization process (recall figure 11). Some remarks are 
given below to highlight the assumed kinetic energy exchange process. 

Communication between primary and secondaiy flow in the present work is based 
in part on kinetic energy gain by the secondary flow due to kinetic energy loss of the 
secondary flow. Such calculations are for the specific purpose of updating the secondary 
flow as time evolves (changes in the primary flow are effected through solution of the 
traditional momentum and energy equations for the same). Three steps are involved: 

1. For a flow condition given after the primary flow has advanced forward in 
time, compute the jet streamline position at station N, for simplicity this 
station is taken to be the point of mixing region discharge. 

2. Compute the entrained flow kinetic energy gain before and after the change 
in primary flow conditions; the difference is the net change in entrained flow 
kinetic energy due to changes in primary flow. 

3. Update the entrained influx kinetic energy for the net change within the 
mixing region; extract the new secondary flow inlet velocity. 
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7. DISCUSSION OF RESULTS 


Characteristic Test Case 

In the absence of data from transient flow ejector experiments (or even from modern 
multi-dimensional Navier-Stokes solvers), "verification" that the proposed ejector analysis can 
provide reasonable thrust predictions must be left to engineering judgement. Because of this, 
a "familiar" ejector forcing function must be used. In the present work the system response to 
a step-function input is not only a characteristic transient case study, but the scenario is also 
coincident with typical STOVL ejector application. 

For demonstration purposes the ejector system response to a step-change in primary 
nozzle flowrate from 18.7 to 21.85 lb^/sec is chosen because (a) experimental steady-state 
data at each of these operating points is available, and (b) the 17% change in primaiy flowrate 
is well beyond a "smalr'-perturbation examination (this exercises the system non-linearities). 
Changes in primary flow stagnation temperature are taken from data given with flowrate data, 
but the corresponding static pressure is computed by the procedure given in Appendix J. 

Calculation Results 

For the mixing region finite-volume length of 0.18 ft and a characteristic mixing region 
velocity of 500 ft/s, the characteristic time step for computations is 

0-18 

A t = s0.4ms 

500 

To avoid infringing on this stability limit a computational time step of 0.1 ms was established; 
100 time steps provided the necessary interval for examination of the step-function test case. 

The empirical coefficient in the transient analysis, Cj (required for calibration of the 
primary-to-secondary kinetic energy exchange mechanism), was selected to match the asyp- 
totic transient thrust prediction with the quasi-steady value at the new set point; satisfying this 
condition required 

C) = 0.97 

Figure 20 illustrates the predicted ejector thrust profile for this coefficient. It appears the 2 
millisecond residence time of the flow (elapsed time for primary nozzle flow to reach the 
diffuser exit plane) is slightly less than the 3 millisecond interval for the thrust to reach a new 
maximum. Oscillations in thrust after that point appear to settle in about 5 
milliseconds. 

An unexpected feature of the thrust profile is the dip in thrust immediately following 
the step-change in primary nozzle efflux. Examination of the field variable profiles reveals 
this is not a numerical problem, but that the increase in static pressure associated with the 
instantaneous change in driving flow temporarily impedes the secondary flow. After a short 
period, the secondary flow kinetic energy builds (commesurate with the increase in primary 
flow energy) to overcome this effect, then continues in the intuitively expected manner. 
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Remarks 


A distinctive second-order flavor is displayed by the predicted thrust profile; under a 
second-order assumption the ejector test case has approximately a 0.75 damping ratio and a 
natural frequency on the order of 300 Hz. Although tne results seem reasonable, it is necessary 
to conduct more extensive computational tests before conclusions about the order or linearity 
(about the perturbation) of the system can be made. The purpose here is to establish the 
routine is operational and that it can in fact provide reasonable results between two steady- 
state conditions. 

The strongest criticizim of the proposed method of analysis probably lies within the 
entrained flow prediction by kinetic energy exchange; the one-dimensional flow limitation 
has required the traditional theoretical analysis of the problem to be modified and an empirical 
coefficient introduced. This may mitigate the robust nature of the simulation approach and 
require fine-tuning for a specific ejector configuration. Once this has been established, 
however, the simulation permits characterization of the transient ejector behavior very quickly, 
for a wide spectrum of operating conditions. If it determined the ejector time constants are 
"small enough" so as to be neglected (relative to time constants for other propulsion system 
components) then a quasi-steady flow assumption for the ejector mixing-region may be valid. 
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8. CONCLUSION 


Assumption Highlights 

It should remain clear the intent of the proposed finite volume method of analysis 
is to meet the combined requirements of 

1. Thrust prediction for real-time simulation, 

2. A predictive fluid-dynamic methodol ogy, and 

3. Characterization of turbulent flow. 

As a result of the assumptions and compromises that must be made to reach these goals, 
two empirical constants nave been introduced. The result is a simulation detailed enough 
to allow a rational introduction of experimental data in the simulation, while at the same 
time being of a simple nature; this is anticipated to provide a realistic candidate for 
real-time simulation. Several assumptions made with this goal in mind (and with significant 
impact on the algorithm structure) are worth repeating: 

1. Primary nozzle and all free-stream conditions are known as a function of time. 

2. Quasi-steady flow conditions exist at the inlet and diffuser. 

3. Entrained flow velocity predictions are adequately given by the proposed kinetic 

energy exchange mechanism. 

Exploration of a step-function test case reveals that items 2 and 3 appear not to have 
compromised the fundamental description of ejector physics. 

Closing Remarks 

The method of analysis for the description of transient ejector characteristics pro- 
vides reasonable results for the single test case considered. As such, the method can be 
concluded as viable if the specific intent of development is kept in mind. An operational 
computer program has been based on the equations presented. Application of the method 
to a broad class of ejector configurations and operating parameters will require exploration 
into the sensitivity of the two empirical coefficients introduced in the course of the analysis. 

As more detailed experimental and theoretical treatment on the theory of turblent 
mixing evolve, it remains to make practical application of those results to the ejector 
analysis. It remains that the basic turbulent flow control volume approach is an excellent 
vehicle on which to test new (or modified versions of the present) method of turbulent 
mixing characterization. 

Solutions to some of the equations described in this work have, due to the limited 
time available, been solved iteratively. There is considerable opportunity to decrease the 
execution time of the method by replacing the iterative rotuines with analytic solutions. 
Such solutions may be available in future work through application of, for instance, the 
MACSYMA symbolic manipulator. 
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Figure 2. Generic Vehicle Flight Control Architecture and Controls 





Figure 3. Typical Ejector Wall Pressure Distribution 
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Figure 7. Nomenclature for Arbitrary Control Volume 
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Figure 8. Flow Model for High and Low Secondary Flow Rates 
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Figure 9. Transverse Velocity Distribution 
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Figure 10 Control Volume Nomenclature for the Ejector 
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Figure 12. Algorithm to Accommodate an Ideal Gas Assumption 
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Figure 14. Nomenclature for Finite Volume Analysis 
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Figure 16. Jet Flow Pattern Approximation 
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Figure 19. Unsteady Flow Solution Procedure 
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Figure 20. Results from Transient Flow Test Case 
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Appendix A 

Development of the Integral Equations 


Overview 

A material volume of V(t) and surface S(t) can be described by four laws of motion: 

1. Conservation of Mass: Based on continuity conditions, this is a statement that the 
material volume is of constant mass. 

2. Balance of Linear Momentum: Newton’s second law provides that the rate of change 
of momentum is equal to the sum of body and surface forces. 

3. Balance of Energy: The first law of thermodynamics describes the balance of internal 
plus kinetic energy, the rate at which work is being performed, and the rate at which 
heat is transferred. 

4. Creation of Entropy: The second law of thermodynamics dictates which of the energy 
transport processes (that the first law provides) that are acceptable, the rate of 
entropy creation is balanced by the sum of the rate of increase of entropy and that 
transported through the material surface. 

These laws are not derivable from a set of more primitive laws and are essentially axioms 
supported by experimental work. 

In the sections to follow the development of the control volume equations is given 
for an arbitrary control volume with a uniform (but time dependent) motion in space. 

Continuity 

Conservation of mass for a material volume takes the form 


— P dV = o 
dtJv 

where d/dt is the total derivative. From Reynolds’ transport theorem 

f ^ 


(/l.l) 


IX" 


dV 


dt 


d V 


p u • /.i d A 


(A. 2) 


Reynolds’ theorem can also be applied to a control volume defined by a region D and 
surface f, 


d_ 

dt 




+ 



U • nd A 


(/1.3) 
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where ^/defines the surface velocity of the control volume, related to the absolute velocity 
by 

V = « - U_ {A A) 

Except for the total derivative terms, integrals over the region D are identical to 
those chosen to be instantaneously coincidental with the (arbitrary) region V; in other 
words 


f^dv - f d ±av 

J v dt J n dt 


( A .5) 

dr d r 

— pdV # — pdV 
dl J v dt J n 


(.4.6) 

but through substitution 



( pdV = ( pdV - <£ pu- nd A + 

dl J n dt J v J v ~ ~ 

£ pU ■ ndA 

(A 1.7) 

so now 



f pdV + <f p(u- U) • ndA = 0 
dt J n J r — — — 


(A 1.8) 


Momentum 


The net body force and surface forces acting on the material volume is balanced by 
the rate of change of the material volume momentum 


d 

dt 


L 


pndV 



p gdV 


n • SdA 


{A. 8 ) 


here, g is identified as the specific gravitational force; the general fluid stress tensor, «§, 
can be decomposed into a "pressure stress" and a "viscous stress" (deviatoric stress) tensor 

5 = ~p[ + S' (A.9) 

To obtain the momentum balance for a control volume, Reynolds’ transport theorem is 
first applied to the total derivative for a material volume 

^lf y pi - dV = X^7 (p - )dl/ + (/i. io) 

and also to the control volume 
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d 

dt 



udV 


' y t (pu)dV 
not 



n)(lA 


(A. 11 ) 


where, again, we note the use of U in the latter. Once the arbitrary material volume is 
extended to coincide with the control volume 


d 

dt 


j" p a d V + j> pu(u • n)dLA - <j> p u((J ■ n)d A ( A . 1 2 ) 


— f p udV = 
dtJv ~ 

It is now evident that 

-j- pudV = -j>pu(u-U)- ndA + J^pgrZF + - d/t (AA3) 


and therefore 

■ f pudV = - j> pu(u-U)- ndA 

! J n 


d_ 
dt , 


f pv 

J 


£'-i 


gdV 


( pl)dA + 




dA {A. 14) 


Energy 

Because textbooks and the literature are often inconsistent in the use of the term 
"energy equation", some clarification is worthwhile so that the present work is interpreted 
correctly Obviously, an energy-conversion law must account for all forms of energy within 
and across the control volume. Confusion often arises since a continuous flow of mass 
across a control volume admits two independent energy equations. In one case, the 
mechanical power equation is derived from the product of the Navier-Stokes equation and 
the oarticle velocity; this gives the so-called balance of mechanical power or transport 
equation for kinetic energy. The other case involves the first law of thermodynamics and 
is the general power equation, or, as viewed here, the general energy equation. It can be 
shown that the heat equation is the difference of the mechanical and general power 
equations; in the present work the heat equation can be transformed into a statement o 
the entropy balance and will be covered in Section A. 5. 


With the first law of thermodynamics in mind, the energy per unit volume of a fluid 
is the sum of the specific internal energy, e, and the specific kinetic energy, u 2 /Z For a 
material volume this energy is balanced by the energy due to heat and work; the genera 
power equation is therefore 




p g • u d V 


+ 


<|> u • ( n * S )dA 

V S “ “ — 


£ 


q -nd A 


(.A. 15) 
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Employing the logic used earlier in the extension of the material volume to a control 
volume 


d 

dt 


fA e O di/ ■ jtiA e A) dv * 


(/1 . 16) 


so that 


d f 

dtJv 


p\ e + ^u.*' \dV = | pg-udl/ + 


ndA 


u • (n • S)dA - £ q • n 

- j> p(^e + ^u 2 ^(u- ndA 


From the definition of the stress tensor we expand the stress term: 

u ■ (n- S)d A = £ u - (n- S')d A - j> pn ■ ud A 

and since the enthalpy is defined as 


( d . 1 7) 


( 71 . 18 ) 


h = q + — 
P 


(71.19) 


then 
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u 
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P 
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{a - U)- ndA 


p g ■ udV 



pn ■ ad A 


£u-(n- S')dA 
j> q ■ ndA 


( A .20) 


It is quite useful to assume the heat equivalent to the work of viscous forces to be 
immediately transferred to the region where the forces occur. No part of that energy is 
transmitted to the surroundings so its presence becomes nested within the other forms of 
energy the governing equations provide a balance for. Furthermore, if gravitational effects 
are neglected and like terms in the energy equation are eliminated, the energy equation 
becomes 
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57/„ p (' l 4-p) d ' / = -{pi- n - dA - 

(A .21) 

Alternate Energy Equation 

In the absence of dissipative effects, an adiabatic ejector representation of the heat 
equation is 

— f podV = - <£ pit ■ ndA (A .22) 

dtJv Js ~ ~ 

Apply Reynold’s transport theorem to the total derivative for a material volume 

— f (pc)dl / = f —(pp)dV + <f> po(u- n)dA (/1.23) 

dtJv Jvot j .*>• 

and also to the control volume 

— f (pc)dl/ = f — {pe)dV + *t> po(u ■ n)dA (A. 24) 

dt J n Jndl .'i 

then 

— f (po)dV = f (pc)dV + 'f po(u ~ £/) • ndA (A. 25) 

dtJv dt J n *s 

and 

— f (pp)dV = - <£ po(u-U)- ndA + f p(un)dA (A.26) 

dt J a Js - — - j s 

Since the enthalpy is defined as 

h = g + l - (/I -27) 

then by rearrangement and substitution 
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p ^ /i - — J ( u - U ) • n d A 


' phv ■ n - f p{/ • ndA 

s ~ ~ J s ~ ~ 


pu - n d A 


( A .28) 


An interesting observation extends from the steady-flow zero-reference velocity 
simplification; the result is 


0 = - phu ■ ndA 


(A. 29) 


which is simply a statement of conservation of heat content between two streams that mix. 
Compare equations A.28 and A.21. 


Introduction of the ideal gas approximation into the heat equation yields 


JtL P { CpT ~ f) dV = 


pU • ndA 


(A. 30) 


Now, from the relations 

p = p RT 


7T f ( Pl dl/ 

dtJnV y- 1 J 


which simplifies to 


' p -v ndA - <p pU ■ ndA (/1 .31) 

s y- 1 - - Js - - J 


dlJ n l Y ~ 1 


d V = 


p -v-ndA - 4> pU • ndA 

Y “ I “ ~ Js ~ ~ 


(3.32) 


and if the specific heat ratio is constant 


■j^pdV = -j> pyv- ndA - (y~\)j> pU_- ndA 


(3.33) 


In the finite volume analysis of the present work this form of the heat equation is more 
convenient (than the general power equation) because of the absence of the cube of 
velocity. 
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Entropy 

Creation of entropy in a material volume is given by the balance of the rate of increase 
of entropy and the entropy flux across the surface and, by definition, is a positive quantity, 


— psdV + 

A • 




d/1 > 0 


( A .34) 


where the entropy flux has been identified as the quotient of the heat flux and temperature 
of transport (this can be shown; see for instance, Appendix B). Reynolds’ transport theorem 
is again drawn on to provide the extention of a material volume to a control volume with 
local surface velocity U.', the result is 


(.4.35) 


-j- J psdV + j> p s {u- U_) • nd A ~ y —q- ndA > 0 (.4.35) 

An interesting analytic excursion (see Appendix B) shows the production of entropy is 
given by viscous dissipation and thermal effects 


MI 


—pV 
p I 


ycf V/' 


PV T 


where the viscous dissipation is defined as 


= (S' • V) • u 


(A.37) 


It is evident that if the entropy balance is non-zero the result must be positive since the 
production terms are proportional to the square of temperature and velocity gradients. 
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Appendix B 

Derivation of the Heat Equation 


Overview 

It is understood that the mechanical power equation is given by the product of the 
Navier Stokes equation and velocity, and that the general power equation is formed by 
application of the first law of thermodynamics to a control volume. Here, an integral form 
of the heat equation is derived as the difference between the integral forms of the 
mechanical power and general power equations. 

Mechanical Power Equation 

For an infinitesimal fluid volume the product of velocity and the Navier Stokes 
equation yields 

du 

PU • -77 = P U-g + W-(V-S) (fl.l) 


where the general stress tensor is composed of normal (pressure) and viscous components 

S = -p/_ + S' (5.2) 

If gravitational effects are ignored and the total derivative written in explicit form, the 
mechanical power equation becomes 

du 

pu- ~77 + P 1 {(« - V)u)} = - ( u • V ) p + u-(V-S) (5.3) 


Expanding the first term 

du _ 3 f pit 2 \ u 2 du 

p - dt ~ J ~?.!n 

then the convective flux term becomes 

pu • {(« • V)u) = (“ ,v )(pj)-j(« , v )P 


(5.4) 


(5.5) 


and combining with the mechanical power equation results in 




this can be reduced to 

= -V-|u^p^j}-(u- V)p + u-(V-S') (510) 

which, when integrated over a material volume V, yields 

Ui{4) av ■ 

(5.11) 

This is the desired result; it is a statement that the integrated kinetic energy inside the 
material volume will change if work has been done (due to pressure forces or viscous 
stresses within the control volume) or there is a net flux of kinetic energy. 

General Power Equation 

The first law of thermodynamics for a system has been shown to be 

p^e + %^j dV = f ^ (n- S)dA- j> q- ndA (5.12) 

If the gravitational effect is ignored and Reynolds’ transport theorem applied to the total 
derivative, 
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[- 

Jvdt 


p\ e + ^\)dV - 


■ j> P^c + l ~~ ’ «d/l + j’ n ■ (« • S)dA - j> q ■ ad A 


(£. 13 ) 


which can conveniently be written 


/i'T 


dl/ 


at 


(pp)dl/ = - ^ ^ jn • ud /1 - ^ pen • ucf ,4 

-fv-qc 

J v ~ 


+ V-(t/- 6 ')dl/- V-gd /1 


(£• 14 ) 


Now, since 

V-fa-S) = u-(V-S' - a • V /’ + <t>-P(n-V)u 

then 


( 5 . 15 ) 


XX( p l)' il/ * f ^ (pc)d|/ ■ 


•at 


/X 


, p — /x • (x d /I 

s V 2 - 


- p o a • 


ad A 


+ (a • (V • S' ) + <!>- xx- VP}dI/ 


V- qdb' 


- f Pn 

J v ~ 


(V- a)dl/ 


(£. 16 ) 

This is the desired form; it is evident we can identify the mechanical power equation within 
the general power equation. 


Heat Equation 


Recall that the heat equation is formed as the difference of the mechanical power 
and general power equations; for the form of these equations as derived above, the result 
is 


at 


(pe)dl / = -<j>ppn-udA + J <| )dV 


y-qdV - J' pn • ( V • a)dV {BAY) 
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or simply 


4- f p edV = f {$-V-q-PS7-u}dV 
dtJv J v ~ ~ 

An important form of Reynolds’ transport theorem provides that 

— pXcLV = f p—dV 
dt Jv Dt 

and therefore 

f^ + v.Z-i + rvujdv = o 

so for an arbitrary material volume 


p— + PV-u = d> - V • f/ 
¥ Dt ~ - 


Introduce 

ID p 

v Dt 


- V • u 


SO 


r 2l = ± + p ! 2l 
Dt Dt Dt 


p { T^-P°> 
' Dt Dt 


Dp 

+ p £ — > 

Dt 


... Ds 

= p / — 

Dt 

= <j> V • q 


Recognize that 

q = - k V T 

and 

kVjT 

T" 


K I 

— + -V-(kVV') 

V 2 /' 


(«.I8) 


(D.19) 


( £ . 20 ) 


(D.21) 


(£. 22 ) 
(£.23) 


(£.24) 


(£.25) 


(£.26) 
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Appendix C 

Steady-State Nozzle Coefficient 


Flow coefficient 

At an (isolated) ejector nozzle exit the steady-state mass flowrate is, for a uniform 
flow, given by 

rh = pu A (C . 1 ) 


Of interest here is a more convenient form of this relation which introduces stagnation 
field variables and the flow coefficient, 4>. 


Stagnation conditions are incorporated by the product 


rh 


p 0 c 0 A 


P_u c\ 
PoCCo/ 


(C.2) 


and identify the flow coefficient 

rh = p 0 c 0 A 4> 


(C.3) 


The flow coefficient can be expressed in terms of the Mach number 

4> = M— — (C.4) 

P o c o 

but a more convenient relation uses 

r ll 

M — = — (C.5) 

Co c 0 

so the expression for the flow coefficient becomes 

1 



To complete the desired form of the flow coefficient we recognize the isentropic pressure 
ratio 


P_i 

Po 



(C.7) 
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and therefore 



As a final step, introduce 


1 2 



(C.8) 


(C.9) 


(C.10) 


Aerodynamic choking 

Channel flow is considered aerodynamically choked if the normalized mass flow rate 
through a section reached a maximum. Because the fluid stagnation properties are assumed 
constant, the flow function has the functional form 


= 



(C.l 1) 


(and represents the functional form of the normalized mass flow rate). 

If is a maximum so must also its square; we obtain the condition for choking that 


d(4> 2 ) 

rf(P 2 ) 


(C .12) 


where 


P = 



(C.l 3) 


The square of the flow function is now 

* 2 - 


(C .14) 
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which simplifies to 


4> 2 


2 

Y- 1 



2 

Y" 1 


P*" 


and the derivative is 


d(<t> 2 ) 

d|3 


0 


2 — 
— p Y '' 
Y- r 


2 

Y“ 1 



Since 


( 3/0 
then it is clear 


(C. 15) 


(C. 16) 


(C. 17) 


111 

Y- 1 


P 


2 

Y- 1 


and therefore the condition of aerodynamically choked flow is that 


(C. 18) 


P 


2 

Y + 1 


(C.19) 


or 



(C .20) 


(C. 21) 


When this result is compared with the expression for the local Mach number we find the 
trivial result that sonic conditions are reached when the flow is choked. 


Summary 

Steady-state mass flowrate through the primary nozzle of the ejector is given by 

m, = p 0 c 0 A 1 4> 


where 
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Appendix D 

Alternate Mass Flowrate Equations 


Summary of Relations 

The literature reflects the use of several different forms of the mass flowrate equation; 
the most common are gathered below: 


m = pc A 


Y" l / 2 Y 


p NV-1/V 


(D.l) 


rh = p A 


2y RTJ ] _( P_ 


Y- I / y\ ' 


(D.2) 


AP 0 ( P 




2y 

R(y- i ) 


( 0 . 3 ) 


For the present work, the mass flowrate equation is that given by the derivation in Appendix 
G: 


rh = -£JLf 6 (y,M) 

v R 7 o 
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Appendix E 

Change in Entropy due to Mixing 


Overview 

In situations where multiple solutions for the Mach numbercan be extracted from 
the momentum equation, the "correct" answer is quite often given by that whose circum- 
stances are in collaboration with the second law of thermodynamics. For a steady-state 
flow, the form of the second law takes on a convenient form and is derived below for 
interest and the completeness of this report. 

Derivation 


For an ideal gas, Gibbs’ relation gives that 


dS 


C'f-* 


dP 

P 


(FA) 


If the gas is thermally and calorically perfect, then 

^ ... Y 


R - const 


Y- 1 


and integration of the previous result gives 


(F.2) 


S.-Sp 

R 


T 

I n I — I - 1 1 1 
• 0 


Po 


(FAS) 


Consider steady-state, steady-flow as an example; the integral entropy balance yields 

J p s u. ■ nd A > 0 (FA) 


so from the definition 
d in 


dA 


p u • n 


(F.5) 


then 


-fr 


S dm - 0 


(F.6) 


For one-dimensional flow 
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-S P] rh Pl -S sl m sl + S z rh 2 > 0 


(E.7) 


we obtain 

AS 


fil s 1 ^71 2 v ~ 

= -S P |-S SI — + s 2 — > o 

nip | am pi 


An p i 

Because 

rh sl = ni 2 ~ m r] 

then 

AS „ f ^2 


An p i 

and so 

AS 


Anp, 


\ fn? 

Sp, -S sl l — - 1 + s 2 — > 0 

P1 511 AAl p i ; Anp, 


iti o 

- -(S,,-S sl ) - -1-(S 2 -S sl ) > 0 


A71 


p 1 


By application of the ideal gas relation 


AS 


Rrh 


p i 


In 


Y-l \T 


T 


p i 


s i 


In 


p i 


s i 


In 


T- 


nipi\\y~ i V 7 


-in 


si 


SI 


which simplifies to the form 


AS 


•In 


T e 


SI 


m- 



I'? \ < n -. 


T s i 


AAl Pl 


In 


Rm P[ Y"1 \ 7 r\ J Y 1 V^pi 

This can be re-arranged in terms of the mass entrainment ratio 


pi 


AS 

/? AAl p ] 

AS 


Y l 7 
— In 


Y 


R rn 


■In 


pi 


7 pi 

T pi 


51 ' + — In 


s \ 


Y- 1 
Y 

Y - 1 V m 


Ll) + JL- ( 

T s I ; Y-IV/»ai 



7 '- 


ATI ' 


m 


SI 


■ p I 



1 ' <: 


S l 




/Pi 

111 , 

nip i V /pi 


m 


S 1 


/n 


p i 


Introducing the entrainment ratio 


(F.8) 


(£*. 9 ) 


(£. 10 ) 


(£. 11 ) 



(£. 12 ) 
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( . 1 3 ) 


t* 


in si 
in pi 


then we have the desired result 


Rmri Y-l V r p i 


HY 


In 



-d + H)ln 



(£. 14 ) 
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Appendix F 

Mixed Flow Fluid Properties 




Basic Relations 

Fluid properties for the mixed flow can be estimated through application of Dalton’s 
law of partial pressures to an ideal mixing process. We summarize for convenience here 
relations extracted from the works of, for instance, Addy et.al.[1981] or Minardi[1982]. 

The specific heat is given by 


C 


p. MR 



Cp.sl 
+ 1^ 


c 


p.p 


where the specific heat is related to the gas constant by 


(F.l) 


(Y- 1) 


R 


The equivalent molecular weight is given by a similar relation 


M 


MR 


= M P ( 1 +|i)< 1 + P 


At, 

at 


The specific heat ratio is 


M ,> 


Y MR 


Y p 


Y s (Y/.-l )Vr 

+ y7(y s -DWs 


(/• . 2 ) 


(F. 3) 


(F.4) 
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Appendix G 
Flow Functions 


Overview 

Manipulation of equations is simplified (and the potential for errors reduced) when 
commonly occuring expressions are introduced as functions. In this appendix five 
expressions with physical significance are derived; they are: 

1. Steady flow stagnation temperature ratio 

2. Steady flow isentropic pressure ratio 

3. Steady flow isentropic velocity ratio 

4. Steady state isentropic area ratio 

5. Mass flow function 

Steady flow stagnation temperature ratio 

For a steady inviscid flow with negligible body forces the energy equation for a 
streamline is 


o 

'■w 

II 

+ 

(G.l) 

where hQ is the stagnation enthalpy. This equation does not imply the flow is isentropic, 
only that there is negligible heat transfer across the control surface. Introduce a perfect 
gas so that 

h = c T + const 

(0.2) 

Pv = RT 

(0.3) 

where P is the thermodynamic (static) pressure; also 


c p = y/?/(y-l) 

(0.4) 

c 2 = (dP/d p) s = y {dP/dv) r = y RT 

(0.5) 


where Cp is the specific heat at constant pressure and c is the speed of sound. These results 
provide^ 


u 


c n r + = r /- 7 'o 


{ 0 . 6 ) 
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(0.7) 


C p / 1 U. C p 1 o 

c 2 2c 2 c 2 


so that 

To 

T 


1 u 2 c 2 

2 c 2 c,,T 


l 4 Al2 ^ tY_,) 7 




In functional form this is written as 


(G.8) 


Y - M y.m) ( C 9 > 

Note that no specific assumption of isentropic flow was made. 

Steady flow isentropic pressure ratio 

In a steady isentropic flow there exists a convenient relation between temperature 
and pressure 


To 

f — ) 

(£• 10 ) 

T 

\p ) 


and thus 

V 


To 

( 1 + Y ~ 1 M 2 \ 

(£• 11 ) 

P 

\ 2 / 

SO 



II 

V 

{ / i ( Y • M ) > * ‘ = 

(£. 12 ) 


Steady flow isentropic velocity ratio 

It is useful in the analysis of choked flow to invert the isentropic pressure ratio 
equation so that 
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( C .13) 


this can be rearranged to yield 


u 

c 



(G.14) 


A characteristic pressure ratio is that given by the quotient of the stagnation pressure 
and the static pressure, 


NPR 



( GAS ) 


and we recognize that, as a general rule, P < Pq, then NPR > 1. The isentropic velocity 
ratio is now 


a 

c 


(N PR) y ~ W2y 



= f,(y,NPR) 


(£. 16 ) 


Since the speed of sound can be written as a function of the nozzle pressure ratio, 


V- I 



(£. 17 ) 


then an alternate expression is 


a 

/ 2 / 

— \ 

Co 

V Y- 1 1 

1 -(NPR) 1 '/ 


(£.18) 


Steady-state isentropic area ratio 

For a steady-state flow the condition of conservation of mass yields the mass flux 

ratio 


A _ P . u * 

A‘ pu 

where an asterisk denotes sonic conditions. From this ratio 


(£.19) 
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I 


p p.C. lPo c o 

p u. p 0 c 0 iVI pc 



= / 5 ( V'M) 

Mass flow function 

From the definition of the mass flow rate 

vAP 

m = pvA = — — 



ih .. I V l o £_ 
PA M v R'I'o I c o 


(C .20) 


(G.21) 


(G .22) 


(G .23) 
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and the temperature ratio 



ft (Y.M) (G. 26) 


Summary 

Five dimensionless functions have been derived, each associated with a specific 
nondimensional physical meaning derived above; they are: 
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Appendix H 

Integrals of Self-Similar Profiles 


Introduction 

Specific assumptions for the dimensionless form of the self-similar temperature and 
velocity profiles will avoid cumbersum numerical integration schemes in the final analysis 
(and in the computer code). This appendix provides an assumed form of the various profiles 
for the present work and the results of integrations of them. Some of the more difficult 
integrals have been explored with the MACSYMA symbolic manipulator on the 
NASA/Lewis Vax. 

If we consider the jet expansion characteristics of Figure H.l, then it is evident two 
expressions should be constructed for each self-similar profile, one for Region I and the 
other for Region II. The basic form of the non-dimensional profiles are given by: 

f : 0 < r, < 1 , ?,>l) 

<t>(%) = ( 0 : 1 < < i , i > 1 > (//••) 

l(l-? l,s ) 2 ; 0<^<t f, < l / 

( l: o , t>n 

A(U = < 0 : ?,*<?,<?,, 1 > 

v i ; o < $ < t . | < i ) 

Application to integrals incorporating these profiles are given in the sections that follow. 
For the case of analysis in the present work it is assumed that 


( // - 3 ) 


Basic f; Integrals 

Several of the velocity integrals are common to many of the self-similar profile 
integrals, so it simplifies the presentation to summarize them. First, set 

/, = ( 1 ~ ' ' 5 ) 2 <"•*) 


from which 

f/.cl* = +0.2l)^ 4 - 0.8^ 5 + ^ (W.S) 

= + O.Bf, ’ :i (IJ .6) 


< 07 :• 



(H.7) 


f 


2 1 10^ - 560?; ° + 1 155^- 1232£ + 7/Of, 

f < d ? ‘ 770 


/ 


2 1 l^ 7 -S6^ 55 + 77^ 4 

(!-./,) d$ = — ^ — 


/ 


220^ - 1 120^ 5 5 + 1925^- 1232^ 
1 540 


(/,-/!)<*$ = 


2.5 


(W.8) 

(f/.9) 


E{ Integrals 

Integral evaluation distinguishes between Region I and Region II values, difference- 
residing within the limits of integration in each case. Figure H.l illustrates the variables 
used. 


Region II coefficients 

/• i = r Ad?, = [‘corn + f\o)ctz = i 

j'o ./0 ^1 

E 2 = f\l-A)dt > = f'(0)dH + f^Ud* = i-1 

y o ./ o 

Case where b > B 

F, = f*(l)dS = | 

J o 

= f\o)dH = 0 


(//. 10 ) 

(//.ll) 


( // - 12 ) 


(H . 13 ) 


Region I coefficients 

In region I b is always less than B, but account of the potential core region requires 
a modification of the lower limit of integration (see figure H.l) 


I'-?. 


f‘ ArtS = 

*'n 0 


Ho 


f ' ( I - A ) <7 ?, = f|- 1 
•'n 0 


(//.14) 


( // - 15 ) 
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Fj Integrals 

Region II coefficients 

F x = / 0 * A4>C ^ = //0(/i)^+j\o)c^ 

= + 0.25( 1 ) 4 - 0.8( 1 ) 2 5 + 1 
= 0.45 

F 2 = 1 - A)<i>d£ = f\o)(/ l )df ; * f*(0)(0)dt 

= 0 

Hereafter, elimination of the integrals of ’0’ can be made by inspection. 

F 3 = f*A(l-<P)d$ = ( 1 - / ! ) ^ 

= -0.25( 1 ) 4 + 0.8( l ) 2 5 

= 1 - F [ = 0.55 

F 4 = f\ 1 - A)( 1 -4>)d^ = J q (1 )dt, = |-1 


(W .16) 
(H .17) 

(H . 18 ) 
(H. 19) 


Region I coefficients 

F , = f°( l )( 1 )dr] 

*' ’lo 

= {-ho) + < 1 - °- 8 ( 1 ) 2 " 4 0.25( I ) 4 } + {0> 

- - r | 0 + 0.45 (// . 20 ) 

F = o 
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( H .22) 


= f'( I )( I - 

J o 

= 0.8(1 ) 5! ' s - 0.25( 1 ) 4 
= 0.85 

f 4 = ^(lKijdf, = ii- l 

G; Integrals 

Region II coefficients 

G, = f A4> z d£ = f A4> z d£ 

1 I 0( 1 ) 7 - 560 ( I ) s - a + I 1 55 ( I ) 4 - I 232( I ) z : ’ + 770 I 

770 

243 _ 486 

770 ~ 1 540 

G ? = f A( 1 -<t>fdt, = f'(l-<t>) 2 dt > 

Jo Jo 

l 1 ( 1 ) / - 56( 1 ) 5 S + 77( 1 ) 4 _ 320 _ 640 

77 ” 77 0 ~ 1540 
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Jo Jo 

,J 220( I y - I 1 20 ( I I 925 ( 1 ) 4 - 1232( 1 ) ? :i \ 
l 1540 / 

6 4 = 0 

6, = j\ 1 - A )( I - <f>)^ = J \ I )dK = (C - I ) 

C 6 = 0 


(//. 23) 


(// .24) 


(H .23) 

414 

1540 

(H .26) 
(11.27) 

(// .28) 

(// .29) 
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Region I coefficients 


( 1 )d t] + <i> A (h] 


(//. 30) 


(l-4>) c/11 = — 


= 2 (<t>- 


414 

1540 


(H .31 ) 


(// .32) 


(/y .33) 


g. 5 = (i)dn = 4 - i 


(// .34) 


G 6 0 


(// .35) 


Hj Integrals 

Presentation of the Hs integrals is given in terms of an arbitrary upper limit to the 
integration; this relates to their use in the kinetic energy analysis of Section 5. 

Region I coefficients 


A( 1 -4>) J d?, = 
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II 
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Figure H.l Non-dimensional Mixing Region Profiles 
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Appendix I 

Remark on System of Units 


Introduction 

Few applications of engineering analyses appear exempt from the need to clarify the 
system of units involved in calculations. Typically, confusion extends from translating back 
and forth between force, mass, and pressure; the usual remark is that a consistent system 
of units must be used in analysis. In the present work the Engineering English system of 
units is used. 

This appendix is intended to summarize the 6 systems of units commonly used in 
engineering analyses so that no confusion will exist over the definition of Engineering 
English system for the present work. The difference in each system of units can be described 
by first identifying the fundamental units in each case, then categorizing each system by 
observing: 

- the magnitude of the fundamental units, 

- the choice of the physical nature of the fundamental units 

- the choice in the number of fundamental units 

The magnitude of the fundamental units originates from the metric and English 
systems. Because metric units are established primarily by international conferences, 
English units are related to them by several convenient conversion factors. Conversion 
factors exist within the different English systems and reflects the varied historical devel- 
opment of the overall system. 

All the systems include length and time as defined quantities - the physical ' nature of 
the fundamental units indicates whether mass, force, or both are defined within the system 
In an absolute (also known as physical) system the mass has a defined fundamental unit 
and force units are derived on the basis of Newton’s second law. An absolute system is 
one for which measurements made in terms of the fundamental units are independent ot 
the location of the measurements. In the gravitational system, however a standard weight 
(standard force) is defined and mass units are derived. Note in either of these systems only 
three defined units, known as primary units, are required to define the system, with all 
remaining quantities derived (secondary units). The engineering system of units in unique 
in that both mass and force are defined. In this case the total number of primary units is 
four and consistency of units is provided by the introduction of a universal constant in 
Newton’s second law. 

From the three general systems described above, six specific systems of units can be 
identified- the absolute MKS, absolute CGS, metric gravitational, English gravitational, 
absolute English, and engineering English In each case force, mass, length, and time are 
related by newtons second law 

F = k(ma ) 
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where k is a constant. A discussion of each system follows. 

Absolute MKS 

Also referred to as the International System (SI), mass, length, and time are defined 
by units of kilograms (kg), meters (M), and seconds ( s ), respectively. The unit of force is 
a Newton (N). A 1 Newton force will give a mass of 1 kg an acceleration of 1 m/s 2 * This 
system is the most popular Metric system and was adopted for international use by the 
Ninth International Congress on Weights and Measures in 1948. 

Absolute CGS 

This absolute system is very similar to the absolute MKS system. The unit of mass 
is defined as the gram (g), the unit of length is the centimeter (cm). Force units are derived 
and given in dynes - no abbreviated symbol - where 1 dyne = 1 g cm/s 2 - 

Metric gravitational 

Since weights are measured by the force of attraction that a given mass experiences, 
a standard weight for the metric gravitational system is defined with the introduction of a 
standard gravitational constant g c . The unit of mass, the kilogram-mass (kgm), is derived. 
Using the meter for units of length and the second for time, the kilogram-force (kgf) is a 
defined standard weight for the attractional force exerted on 1 kg mass by the earth where 
the gravitational constant has a standard value ofg c = 9.80665. Note a total of three defined 
units and one derived. 

English gravitational 

Making use of the foot and second that is common to all English systems of units, 
the English gravitational system differs from the absolute system in the same way the 
metric graviational does. That is, the unit of force is defined as pound-force (lbf) and the 
unit of mass is derived as the slug - no abbreviated name. From Newton’s second law we 
obtain 1 slug = 1 lbf s 2 /ft. When a force of 1 lbf is applied to a mass of 1 slug it will yield 
an acceleration of 1 ft/s 2 . 

Absolute English 

This system is also based on three defined fundamental units; feet (ft), seconds (s), 
and pound-mass (lbm). The units of force are derives as the poundal, where 1 poundal = 
1 lbm.ft/s 2 . 

Engineering English 

This system is based on four fundamental units for length, force, time, and mass. As 
before, the units of length and time are feet and seconds respectively, but the unit of force 
is the lbf and the unit of mass is the lbm. Compatability of units is provided by the intro- 
duction of an English gravitational constant g c into Newton’s law. This is the system used 
in the present work ! 
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F 


a 

m — 

Qc 


where g c = 32.174 lbm ft/lbf s 2 - So for this system, rather than a unit force imparting a 
unit acceleration to a unit mass, a unit force causes a unit mass to accelerate with a value 
equal to the gravitational constant. It should be noted that the constant g c is numerically 
but not dimensionally equal to the standard gravitational acceleration g. Two problems 
arise that often lead to confusion with this system: (a) gc is often thought of as g, implying 
the problem depends on terrestrial graviation when it may not, and (b) the ratio g/g c 
obtained by setting a=g in Newton’s law is often approximated as unity and discarded, 
leaving the appearance that the problem is not a function of gravitational effects where it 
might indeed be an important part. 


Conversion Between Systems 


The present work uses the Engineering English system of units. A summary of the 
"hierarchy" of the systems in given in Table 1.1. Comparison of the MK.S and Engineering 
English units for quanities of key interest are given in Table 1.2. More detailed conversions 
between the various systems are provided in Table 1.3. 


Illustrative Calculations 

Some illustrative calculations may seem trivial, but are quite illustrative in following 
the presence (or absence) ofg c when, for instance, the details of the computer program 
are being traced. 


Gas Constant The gas constant for air at "standard" conditions can be computed from 


R = 

R/M 


where 



R = 

49,720 

/ ft 2 

\ R Lb mole . 

M = 

28.98 l 

' lb m \ 

l mole 1 

Now 



R = 

49,720 

j ft 2 lb m 

28.98 

\ S? R mo le 


mole 


lb, 
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= 1715.66 


ft* 

s ? R 


Specific Heat For calculations in the computer programs the specific heat at constant 
pressure is given by 


so for the data of the previous calculation and a specific heat ratio of 1.4 


c 


p 


1 .4 

1.4-1 


1/15.66 


6004.81 


ft 2 

s 2 R 


but this does not immediately conjure the "0.24" value one might expect. That value is 
obtained by conversion of units 


c 


p 


6004.8 


ft 2 

s 2 R 


I Bln 1 Lb ,s ? 
7/8.6 Tb } ft g~ c lb m ft 


0.24 


Btu 

I b rn R 


Sound Speed For an ideal gas approximation the speed of sound is given by 

c = 

so for the data 
Y = 1.4 

T = 759.67 R 
R= 1715.66 ft 2 / Ra 2 


then 


c 


V 1 .4(759.67) ( 1/18.66) 



= 1 350.80/</s 
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Gas Density Consider an ideal gas at a pressure of 4233.6 lbf/ft 2 and a temperature of 
300 F; the density is determined from 


p = P/RT 

4233.6 / IfryZ/t 2 \ 

1715.66(759.67) \(ft 2 /s 2 R)R ) 

lb ,s 2 

= 0.0032483 — — 

ft' 

To be in a consistent Engineering English system of units the factor of g c must not be 
dropped 


P = ( P/RT)g c 


= 0.0032483(32.174) 


= 0.1045 



Of course, in the English gravitational 


1 slug 


lb yS 2 

7 


l ft 4 lb ,s 2 j 


system the derived mass is the slug 


so the first result is immediately recognized as 


P 


0.0032483 


slugs 

77 


Stagnation Pressure The stagnation pressure for a 300 F gas flow with a 150/0 velocity 
and 4233.6 lbf/ft 2 static pressure can be computed from 


so 


P 


0 


4233.60 


ibA 

ft 2 ) 


0 . 1 045 ( 1 50 ) 2 

4 - 

2 


U> m ft 2 \ 
ft 2 s 2 j 
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P 0 = 4233.6o(^M + 1 175.63/ — 

\ft 2 S \ ft- s' 


and the inconsistency in units resulting from the ommission of g<.; correcting the density, 


= P 


p a 

2 g c 


= 4233.60 + 


= 4270.14 


1 175.63 
32.174 

lb, 

It 2 ' 


Lb 

1 T‘ 


Note that computations in the Gravitational English system do not, by default, require 
an explicit inclusion of a numerical value forg c since it is unity. 


Stagnation Temperature Correction of the density for stagnation temperature calcu- 
lations is done the same as in the stagnation pressure case 


Kinetic Energy of the Flow This quantity has the units of work since the "mass" in this 
case is a flowrate term 


* 

2g c 
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DIMENSIONAL SYSTEMS 



\lb„ - 0.453592fcg 1 lb, - 4.448W 

1 slug - 32.174ft> m 1 ft - 0.3048m 










Table 1.2 Summary of Engineering and MKS Systems 


Quantity 

symbol 

English 

MKS 

Length 

L 

ft 

m 

Time 

t 

s 

s 

Mass 

m 

torn 


|j Force 

F 

lb or Ibf 

N 

Density 

P 

tbm/ft 3 

kg/m 2 

Dynamic 

Viscosity 

P 

Ibf-s/ft 2 

Ns/m 2 

Kinematic 

Viscosity 

V 

ft 2 A 

m 2 /s 

Pressure 

p 

ibf/ft 2 

N/m 2 

Work, Energy 

w 

lbf-ft 

N-m 

Power 

w 

lbf-ft/s 

N-m/s 


Note: 

Qc 


32. 174 


lb f s 2 


kq - m 

1 — 7 

N -s 2 
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Table 1.3 Conversion Factors 
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Appendix J 

Computation of Nozzle Exit Static Pressure 


Introduction 


Instrument data for the primary nozzle flow conditions has been assumed to be for 
mass flowrate, stagnation temperature, and static pressure. From this information the 
velocity, static temperature, and static pressure for the primary nozzle can be computed. 
However, if stagnation pressure is given rather than static pressure the analysis for primary 
nozzle conditions (given in the main body of this report) must be modified; the details for 
this are given below. 

Derivation 


The nozzle pressure ratio can be given by 


Po 

p 


] - 


v~ i 

2 


M ' 


Y/(Y-1) 


so the Mach number is 



Combine this result with the modified mass flow-rate expressions of Appendix D to obtain 


W 

A 


P o 


nmi 

V RT oY- 1 \ 



Y- l/Y 

- 1 




-y 1 /2 Y 


and after re-arrangement 

/vy Y -i/?T 0 _ (Po 
{ A J 2 y P ? 0 V P 



y •* I 

Y 


If the stagnation pressure, stagnation temperature, and mass flowrate are known, then this 
equation contains only the pressure ratio as an unknown; in polynomial form 


P 

P~o 



W V y- 1 RT 0 

AJ 2y ~Pl 
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Solution 


The root of the equation above provides the solution for the nozzle pressure ratio 
and since the stagnation pressure is known, the static pressure computes directly there- 
from. Because this is not an ordinary polynomial which benefits from a simple analytic 
solution a numerical approach is reccommended. Figure J.l provides a listing and sample 
calculation for the APL computer program that provides this service. 
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Figure J.l APL Solution Approach for Static Pressure 
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PO 
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4 
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C53 

TO 
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4 
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4 
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[73 
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i 
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. . . INPUT )A7A. . . 


PO 

= 

2799. 9000 

LBF.- FT2 

TO 

= 

769.7000 

: R 

A 


. 4200 

FT2 

GAMMA 

= 

1 . 4000 


R 

= 

52. 3029 

FT- LBF/ LEM 

MDOT 

= 

IS. 70 DO 

LEM/ 2 
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p 


2219. 7B00 

LET. FT2 
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Appendix K 

Overview of some PLF ejector data 


Some geometric and instrumentation data 

In the evaluation of the methods of analysis outlined in the text, it is necessary to 
have data with which to compare theoretical predictions. Information on some past 
NASA/Lewis PLF ejector tests have been made available to the present work. Figure 1 
summarizes some geometric and instrumentation data believed to apply to the DeHaviland 
ejector tests conducted at the NASA/Lewis PLF facility. 

Knowledge of the type of test data to be collected is fundamental in shaping the 
input requirements of practical computer programs. For the present work the requirements 
developed in the text (see chapter 5 of the text) have been based on the type of information 
one might expect from the PLF. 

The general geometric characteristics of the DeHaviland ejector are also shown in 
figure 1. Use of this data has guided the geometric data requirements for the analysis as 
illustrated in Figure 1 1 of the text. 

Some performance data 

Performance data used for comparison with the theoretical predictions is shown in 
Table 1 - 3. The nomenclature for the computer printouts is defined (to a limited extent) 

in Figure 1. 
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Figure K.l Some Geometric and Thermodynamic Ejector Information 
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Table 1. Data Summary for Preliminary PLF Tests 


Reading 

F-S Temp 

NPR 

Thmst 

PHI 

m | 

223 

769.70 

1.35 

945. 

1.86 

18.70 | 

239 

760.60 

1.46 

1231. 

1.87 1 

21.85 1 

225 

761.70 

1.70 

■a 

1.80 

27.27 

237 

763.50 

1.82 

- 


1.79 

29.66 

227 

763.50 

2.06 

2471. 

1.78 

34.26 

235 

765.50 

2.27 

2910. 

1.78 

37.97 

229 

763.50 

2.38 

3066. 

1.74 

39.96 

1 

233 

758.40 

2.56 

3458. 

1.76 

43.22 

231 

764.80 

2.79 

3839. 

1.72 

46.97 
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Table 2. Data from PLF Test #223 


NASA- LEWIS 


PRELIMINARY DATA 


11/30/87 


PDP11TWO 


REC 07/16/87 19:08:50.656 


FAC PLr 


PGM DO 0 1 


POWERED LIFT FACILITY SINGLE-SIDED EJECTOR FOR SPEY/E7 MODEL RUN 9 


EJECTOR NOZZLE TYPE NOTCHED- CONE NOZZLES 


FORWARD DUCT VALVE ANGLE 90.0 DEG 


OPTION 1 
RUN 

SUMMARY 

FDVA 

NPR 

PRIN 

PR6 

TP 

CYCLES 

9 

90.0 

1.354 

1 .426D 

1 . 370 

769.7 

10 

H 

WNOZ 

CDN 

XM 

XM2 

XIE 

XI6 

19.51 

18.70 

1.0435 

945 .D 

938 . 

532. 

541 . 

PHINET 

PHINET2 

ANET6 

XMSTD 

H 

PA 

TA 

1.7776D 

1.7651 

1.7469D 

96 7 .D 

5.97D 

14.36 

542.8 

WF 

PR 5 

XI5 

ANET5 

XIIN 

ANETIN 


0.05 

1.378D 

546 .D 

1 . 7 308D 

573 .D 

1 .6491D 



EXIT-RAKE 1 LOCATION 
EXIT-RAKE 2 LOCATION 


/L = W. Si 

~P Q = /•?>'•? 7L 


0.00 

0.00 


2 m u 


NASA-LEWIS PRELIMINARY DATA 
POWERED LIFT FACILITY 


11/30/87 PDP11TWO REC 07/16/87 19: 08:50^656 


FAC PLF 


~ • - A * O A LVLU 

EJECTOR NOZZLE TYPE NOTCHED-CONE NOZZLES 
OPTION 4 EJECTOR-NOZZLE CALCULATIONS 


SINGLE-SIDEDEJECTOR FOR SPEY/E7 MODEL RUN 9 

FORWARD DUCT VALVE ANGLE 


90.0 DEG EXIT-RAKE 1 LOCATION 
EXIT-RAKE 2 LOCATION 


PGM DO 0 3 


0 . 00 

0 . 00 


PN 

01 

18.44 

02 

18.82 

03 

18.88 

04 

18 . 76 

05 

18.66 

06 

18.99 

07 

18.98 

08 

18.80 

09 

18.92 

10 

18.69 

11 

18.85 

12 

18.70 

A 

18.79 

PTN 

19.18 

19.27 

19.34 

19.40 

19.44 

19.48 

19.49 

19.54 

19.54 

19.57 

19.56 

19.58 

19.45 

PRTN 

1.336 

1 . 342 

1 . 346 

1.351 

1 . 354 

1 . 357 

1.357 

1.361 

1 . 360 

1.363 

1.362 

1.363 

1.354 

DELPTN 

0.61D 

0.52D 

0.45D 

0 . 39D 

0.35D 

0 . 31 D 

0. 30D 

0.25D 

0.26D 

0.22D 

0 . 23D 

0.21D 

0.34D 

DELPTN/PT5 

0 . 0 3 ID 

0.027D 

0 . 02 3D 

0 . 02 0D 

0 . 01CD 

0 .016D 

0. 015D 

0. 012D 

0 . 01 3D 

0 . 01 ID 

0.012D 

0. 011D 

0 . 017D 

DELPTN/ (PT5- PA) 

0.113D 

0. 097D 

0. 084D 

0 . 071D 

0. 064D 

0. 057D 

0.055D 

0.045D 

0. 047D 

0 . 140D 

1 . 042D 

0.039D 

0 .063D 

DELPTN/Q5 

0.821D 

0.7C5D 

0.61 ID 

0 .521D 

0.466D 

0.416D 

0.400D 

0 .332D 

0 . 344D 

0.293D 

0.310D 

0 .286D 

0.459D 


ORIGINAL PAGE I'j 
OF POOR QUALITY 
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Table 3. Data from PLF Test #225 


OF POOR QUALITY 


NASA-LEWIS PRELIMINARY DATA 


1 1/ 30/87 


PDP1 1TUC) 


REC 07/16/87 1 9 : 1 3 : 08 . 544 


FAC PLF 


POWERED LIFT FACILITY SINGLE-SIDED EJECTOR r0R ^SPEY/F7 ^MODEL 
EJECTOR NOZZLE TYPE NOTCHED-CONE NOZZLE-. 


RUN 9 

rOPWAPr DUCT VALVE ANGLE 90.0 DEG EXIT-RAKE 1 LOCATION 

EXIT-RAKE 2 LOCATION 


OPTION 1 

SUMMARY 






RUN 

9 

FDVA 
90 . 0 

NPR 
1 .697 

PRIN 
1 .S01D 

PR6 

1.715 

TP 

761.7 

CYCLES 

10 

U 

27.45 

WNOZ 
27 .29 

CDN 
1 . 0056 

XM 

1725. D 

XM2 
1719 . 

X IE 
96 7 . 

XI6 
976 . 

PHINET 
1 . 7 8 3-5D 

PHINET2 
1 .7 7 7 3 

ANET6 
1 . 76 7 5D 

XMSTD 
1 765 .D 

H 

6 . 46D 

PA 

1^.36 

TA 

540 . 3 

WF 
0 . 06 

PR5 

1.726D 

XI5 
981 -D 

ANET5 
1 .7583D 

XIIN 
1016 .D 

ANFTIN 
1 .t985D 



PGM 


0 . 0 C 

0 . 0 c 


NASA-LEWIS PRELIMINARY DATA 11/30/87 PDPUTtO REC 07/16/87 1 9 : 1 3 : 08 . 544 FAC PLF *" j1 ‘ 

S r “r?Ki;"»E“»L»“I^E .... K. gn ; .« l [oc.tio. .... 

OPTION 4 EJECTOR-NOZZLE CALCULATIONS 

PN 
PTN 
PRTN 
DELPTN 
DELPTN/PT5 
DELPTN/ < PT5 -PA ) 

DELPTN/Q5 


01 

22.65 

02 

23.26 

03 

23.34 

04 

23.17 

05 

23.02 

06 

25.54 

07 

23.54 

08 

23.26 

09 

23.43 

10 

23.06 

11 

23 . 33 

23 

12 
. 08 

A 

23 . 22 

23.92 

24.15 

24.24 

24.30 

24 .33 

2<*. 48 

24 . 51 

24. 55 

24.53 

24.49 

24 . 54 

24 

.51 

24 . 38 

1.665 

1.681 

1.688 

1.692 

1.694 

1 704 

1 . 706 

1.708 

1.708 

1.705 

1.709 

1 . 

707 

1.697 

0.87D 

0 .64D 

0.54D 

0 .48D 

0.45D 

0.31D 

0.28D 

0.26D 

0 .25D 

0 .29D 

0.24D 

0 

. 27D 

0 . 41D 

0.C35D 

0. 026D 

0 . 022D 

0 . 020D 

0.016D 

0 0 1 2D 

0 . 01 ID 

0 . 010D 

0 . 010D 

0 . 012D 

0 . 010D 

0 . 

01 ID 

0 . 0 16D 

0 . 083D 

0 . 06 ID 

0.052D 

0 . 047D 

0 . 043D 

0 029D 

0 .027D 

0.025D 

0 . 024D 

O .028D 

0 . 02 3D 

0 . 

0 2 6 D 

0 . 039D 

0.738D 

0.544D 

0 .46 3D 

0.413D 

0 .385D 

0 26 1 D 

0 . 236D 

0.219D 

0.215D 

0 .248D 

0 . 206D 

0 . 

2 3 1 D 

0 . 34 7D 


D00 3 


DO 0 3 

0 

0 
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Appendix L 

Jet Boundaiy Streamline 
Equivalent Secondary Mass Flow 

In the exchange of kinetic energy between mixing streams, quantifying the kinetic 
energy gained or lost by a single stream requires definition of a dividing streamline for the 
figure 1 marks the dividing streamline with the jet boundary parameter b* . This is 
different than the boundary defined by point / on the figure, which simply marks the vol- 
ume consumed by turbulent jet expansion. 

The dividing streamline easily derives from the mass flow relation 



F i = 0.45 - ( + 0.255 * - 0.8£j' 5 + £ , ) 
F 3 = 0.55-(-0.2S$f + 0.8£f s ) 

F 4 = 2 ~ 1 


Substitution and cancellation of like terms yields 

C l^ +r ,?^f 5 + C 3^1 + r 4 = 0 

where 


ORIGINAL FkGk. i:S 

OF POOR QUALITY 
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1 


Cl 

C 2 


b_ 

b 


0.25(u m . 2 -y.. 2 ) 

0.80(y a2 - v m 2 ) 


c 3 
C 4 


U m . 2 
p j 5 B — b 0 
P 2 b/ 


-u 


0 . I 


B 

-c, -c 2 -c a -| r _ 1 l y 


The root to this equation provides the location of the dividing streamline, b , at desired 
axial location, station 2. 


Numerical Solution 


The polynomial describing the dividing streamline location does not have an conve- 
nient analytic solution, so numerical methods are used. In the present work the modified 
Newton’s method is used. For this method two derivatives are required: 


/(O = c ,^ 4 + c 2 t; 2 ti -*-c 3 t; + c 4 

/'(O = 4c 1 ^' , + 2.5c 2 t > 15 + c 3 

/"(O = 12c , + 3.75c 2 ^° 5 


The iterative method of solution calls for the following steps 

1) Assume an initial value for the root; here t,c-i; 0 = b 0 /b 

2) u(0 = /ttV/'CO 

3) u '(0 = 1 -/a)/"(£V(/'(0) 2 


4) 6 = -u(^)/u'(t) 

5) If 6 < g then exit otherwise refine the approximation for the root = K + & a °d 
repeat steps 2-4. 
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Figure LI. Jet Boundary Streamline 
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